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ABSTRACT 

A  generally  applicable  method  of  obtaining  kinetic  parameters  from 
temperature-programmed  thermogravimetry  is  presented.  Factors  influencing 
the  selection  of  a  particular  method  for  the  numerous  treatments  reported  in 
the  literature  are  discussed  in  detail.  The  method  of  Friedman  involving  the 
use  of  several  thermograms  at  different  heating  rates  and  determining  Arhennius 
parameters  at  each  percent  conversion  was  chosen.  The  e3q)erimental  procedure 
and  a  method  of  handling  thermogravimetric  analysis  (TGA)data  and  calculations 
by  computer  are  fully  described.  The  application  of  the  treatment  to  some 
specific  pol3mier  degradation  systems  is  reported  in  order  to  illustrate  the 
scope  of  the  method  and  its  potential  usefulness  in  obtaining  information  con¬ 
cerning  complex  degradation  mechanisms.  Poly(tetrafluoroethylene),  an  ali¬ 
phatic,  and  an  aromatic  polyamide  were  the  polymers  selected  for  this  study. 
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SECTION  I 
INTRODUCTION 

A  knowledge  and  understanding  of  polymer  degradation  processes  is  nec¬ 
essary  to  improve  the  performance  of  polymers  in  high  temperature  applications 
and  to  direct  research  towards  more  thermally  stable  systems.  One  of  the 
most  important  methods  of  studying  polymer  degradation  is  to  examine  the 
kinetics  governing  the  breakdown  reactions  from  which  information  about  the 
thermal  stability  and  the  mechanism  of  degradation  can  be  obtained. 

In  degradation  kinetics,  the  rate  of  change  of  some  property  must  be  mea¬ 
sured  as  a  function  of  time,  temperature,  or  conversion;  any  of  several  proper¬ 
ties  can  be  utilized  if  the  property  can  be  correlated  with  the  degradation 
process,  e.g.,  mass,  molecular  wei^t,  a  chemical  group,  or  production  of 
volatile  products.  A  common  method  has  been  to  follow  weight  changes,  A 
thermobalance  can  be  ixsed  to  continuously  record  the  weight  changes.  ]h  the 
past,  isothermal  studies  have  been  most  common,  in  which  the  sample  is  main¬ 
tained  at  a  constant  temperature  while  weight-time  measurements  are  made.  In 
recent  years,  temperature-programmed  methods  have  been  increasingly  used. 
In  this  method,  temperature  is  continuously  raised,  usually  linearly  with  time, 
and  a  thermogram  of  wei^t  versus  temperature  obtained.  In  theory,  the  pro¬ 
grammed  method  should  have  certain  advantages  over  the  isothermal  method. 
An  important  criticism  of  isothermal  methods  is  that  the  temperature  of  the 
sample  cannot  be  raised  instantaneously  to  the  desired  temperatiu'e  so  that 
some  weight  may  be  lost  before  the  degradation  temperature  is  reached.  This 
weight  loss  may  give  rise  to  an  apparent  maximum  in  the  rate  curve  which  masks 
the  true  initial  features.  Initial  rate  characteristics  are  very  important  since 
end-group  and  impurity- induced  reactions  may  show  t;?).  Temperature- 
programmed  methods  should  overcome  this  problem.  It  was  further  visualized 
that  since  weight-temperature  dependencies  were  contained  in  a  single  chart 
the  thermogram  from  a  programmed  experiment  would  yield  the  equivalent 
information  of  a  large  family  of  isothermal  e3q)eriments.  Thus,  sample  uni¬ 
formity  errors  would  be  avoided  and  the  procediire  would  be  less  time  consuming. 
The  thermograms  are  also  useful  for  qualitative  comparisons  of  thermal 
stabilities. 
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In  this  report,  some  of  the  numerous  methods  of  obtaining  kinetic  pa¬ 
rameters  from  programmed  TGA  are  discussed.  The  methods  are  evaluated 
against  the  criteria  that  the  method  should  be  generally  applicable,  should  give 
meaningful  kinetic  parameters,  and  should  shed  light  on  the  mechanism  of 
degradation.  The  selection  of  such  a  method  and  its  application  to  some  pol3mier 
systems  is  also  described. 


SECTION  II 
DISCUSSION 


1.  GENERAL  CONSIDERATIONS 

The  classical  kinetic  expression  which  is  widely  applicable  to  gas-phase 
and  solution  reactions  is  represented  by  Equation  1 

=  kF(C)  (I) 

d  T 

C  =  concentration  of  reactant 

t  =  time 

k  =  rate  constant 

F(C)  =  function  of  C 


hi  classical  kinetics,  F(C)  can  often  be  expressed  as  a  power  function,  c”,  for 
which  n  is  defined  as  the  order  of  reaction.  The  rate  constant  is  temperature 
dependent  and  is  defined  by  the  Arrhenius  equation 


k 


/RT 


(2) 


2 


AFML-TR-68-181 
Part  n 


A  =  pre-ejqponential  factor 
E  =  activation  energy 

Si 

R  =  general  gas  constant 
T  =  absolute  tenaperature 

The  normal  kinetic  approach  is  to  determine  rate  constants  for  a  given  reaction 
at  various  temperatures  using  Equation  1  or  an  integrated  form  of  it  and  thence 
to  calculate  the  parameters  A  and  E  from  Equation  2.  The  activation  energy 
can  often  be  correlated  with  the  breaking  of  specific  chemical  bonds  and  gives 
important  information  concerning  the  mechanism  of  the  reaction  being  studied. 

Polymer  degradation  kinetics  are  normally  studied  in  the  solid  or  melt  phase 
and,  since  a  chemical  reaction  is  occurring,  it  is  assumed  that  a  kinetic  treat¬ 
ment  and  rate  expression  is  applicable.  Owing  to  the  complexity  and  variety  of 
polymer  decomposition  schemes,  it  is  found,  however,  that  the  concentration 
of  polymer  molecules  is  not  eqxiivalent  to  the  concentration  of  reactant  in 
normal  reactions  so  that  the  term  “concentration”  must  be  used  with  extreme 
caution  in  discussing  polymer  degradation.  Two  well  established  types  of  deg¬ 
radation  illustrate  this  point  (Reference  1).  hi  the  first,  degradation  is  initiated 
at  a  chain  end  and  proceeds  by  unzipping  throu^  the  entire  polymer  chain.  The 
sample  loses  one  polymer  molecule  and  its  corresponding  weight  and  volume, 
but  the  concentration  of  polymer  molecules  remains  unchanged,  hi  the  second 
type,  the  elimination  of  side-grov^is  in  the  chain  alters  the  chemical  structure 
of  the  poljuner  but  does  not  necessarily  influence  the  number  or  concentration 
of  polymer  molecules,  hi  polymers,  discussion  in  terms  of  reactive  sites  which 
could  be  chain  ends,  particular  bonds  or  groups  in  the  structure,  the  links 
joining  monomer  units,  etc.,  is  often  more  fridtful  for  elucidation  of  mechanisms. 

It  is  found  in  practice  that  most  degradation  reactions  result  in  loss  of 
weight  owing  to  the  formation  of  smaller,  volatile  species  ^  rearrangements  such 
as  the  coloration  of  poly  (acrylonitrile)  are  notable  exceptions  | .  Therefore,  it 
has  been  found  convenient  to  study  kinetics  in  terms  of  wei^t  loss;  to  this 
extent,  the  approach  is  empirical  since  the  only  proof  of  validity  is  that  exper¬ 
imental  results  appear  to  fit  the  adopted  kinetic  expressions.  In  order  to 
preserve  a  close  resemblance  to  normal  kinetic  procedure,  the  weight  terms 
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should  be  expressed  in  fractional  form.  Further,  since  a  reactant  shoxald  have 
zero  concentration  on  completion  of  reaction,  the  weight  term  should  be  cor¬ 
rected  for  any  residue  weight  remaining  after  degradation.  A  general  rate 
expression  based  on  weight  terms  can  now  be  postulated 


I  dW  _  /  W-Wf  \ 

Wq  -  Wf  d  t  ■  ^  V  Wq-  / 

W  =  instantaneous  weight  of  sample 
=  initial  weight  of  sample 

Wj  =  final  or  residue  weight 


k  F(W) 


(3) 


W  -Wj 

TjT — 7^  =  fractional  weight  remaining 
vV  -  w, 
o  f 


W  -Wj  w^-w 

(1  -  w  -w'"'  "  =  fractional  weight  loss) 

of  of 


Equation  3  is  equivalent  to  that  used  by  Doyle  (Reference  2)  except  that  weight 
terms  are  retained  here  in  preference  to  Doyle’s  active  weight  fraction,  h.  It 
must  be  emphasized  that  the  functionality  F(W)  need  not  be  simple  and  should 
not  be  assumed  in  advance.  The  presumption  of  “order”  type  functionality  should 
be  avoided  unless  it  is  warranted  by  experimental  evidence.  Two  particular  forms 
of  F(W),  which  have  been  observed  in  isothermal  studies,  merit  a  brief  dis¬ 
cussion  at  this  point; 


a.  F(W)  =  (for  simplicity,  assume  =  O) 

'^o 


This  is  a  common  case  which  results  in  apparent  first  order  kinetics  by 
analogy  to  classical  rate  laws.  The  rate  of  weight  loss  is  proportional  to  the 
residual  weight  of  polymer.  Polymers  which  degrade  predominantly  by  unzipping 
show  kinetics  of  this  type,  the  weight  loss  being  a  direct  measure  of  the  amount 
of  unzipping  and  degree  of  degradation.  This,  of  course,  is  a  simplified  view  of 
the  overall  kinetics,  and  many  other  interesting  kinetic  dependencies  can  be 
observed  (References  1  and  3). 


b.  F(W)  =  (  L  -  I  ) 


W 


W, 


(  I  -a  ) 


L  (l\]-L)(L-l  ) 
N 


(4) 
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L  =  smallest  chain  length  which  does  not  evaporate  under  degradative 
conditions 

N  =  initial  chain  length  of  polymer 

a  =  l-e“^  (temperatiire  constant)  =  degree  of  degradation  (bond  scission) 

This  complex  form  of  F(W)  has  been  developed  by  Simha  et  al  as  the  solution  to 
the  random  type  of  degradation  mechanism  (References  4  and  5),  This  mechanism 
consists  of  chain  rupture  taking  place  randomly;  the  act  of  scission  need  not 
result  in  weight  loss  which  is  possible  only  when  scission  occurs  near  a  chain 
end  to  give  a  chain  fragment  of  length  <L.  The  rate  of  weight  loss,  corresponding 
to  a  spectrum  of  products  of  chain  length  from  one  to  L,  is  proportional  to  the 
number  of  chain  ends  in  the  system,  which  increases  at  first  owing  to  the  random 
scission  of  the  chains.  Eventually,  the  average  chain  length  becomes  small 
enough  that  further  increase  in  the  nuanber  of  ends  is  balanced  by  loss  through 
evaporation.  The  rate  of  weight  loss  passes  through  a  maximum  and  subsequently 
falls ;  it  has  been  shown  theoretically  that  the  maximum  should  occur  at  ap¬ 
proximately  26%  conversion.  Pol3/methylene  and  poly  (methyl  acrylate)  are 
examples  of  vinyl  polymers  exhibiting  the  above  features.  Many  condensation 
type  polymers  such  as  cellules ics,  polyamides,  and  polyesters  also  show  random 
characteristics  (Reference  3). 

As  has  been  discussed,  several  functionalities  of  weight  are  readily  ob¬ 
served  in  isothermal  weight-loss  studies.  In  contrast,  primary  data  from 
programmed  TGA  studies  is  much  less  revealing  and  differences  caused  by 
changes  in  weight  functionality  may  be  very  subtle.  The  weight-temperature 
thermogram  and  the  derivative  rate-temperature  curve  (Figure  1)  are  somewhat 
general  in  character,  apart  from  the  obvious  difference  between  a  simple  and 
multistep  process ;  differences  in  shape  are  relative  and  could  be  divined  only 
by  reference  to  other  “known”  curves.  That  both  activation  energy  and  mech¬ 
anism  affect  the  shape  of  the  curve  adds  to  the  difficulties  (Reference  6).  As  is 
shown  in  this  reference,  rate-conversion  curves  are  of  greater  diagnostic 
importance,  although  in  practice,  experimental  fluctuations  from  a  smooth  curve 
mi^t  easily  mask  the  diagnostic  featvtres,  namely  the  conversion  for  the  max¬ 
imum  rate  and  the  initial  slope  of  the  curve.  The  ability  of  kinetic  analysis  to 
afford  not  only  the  Arrhenius  parameters  but  also  supplementary  information 
about  the  mechanism  is  investigated  in  the  following  paragraphs. 
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2.  KINETIC  ANALYSIS  OF  PROGRAMMED  WEIGHT-LOSS  DATA 

Analytical  naethods  can  be  broadly  divided  into  those  using  differential  and 
integral  treatments.  Differential  methods  are  based  on  Equation  3.  In  pro¬ 
grammed  TGA,  k  in  this  equation  becomes  a  variable  since  temperature  and 
therefore  the  rate  constant  are  continuously  changing.  The  expression  can  be 
rearranged  as  shown  for  a  linear  heating  rate,  B, 


The  integral  approach  was  developed  to  use  TGA  weight- loss  data  directly  and 
to  avoid  calculation  of  rates.  The  equation  of  the  thermogram  (References  2 
and  7)  is  obtained  by  integrating  Eqxiation  5 


dW 

:  Wq-W^)F(W) 


The  result  is  given  as 


(6) 


(7) 


in  which  p(E  /RT)  is  a  complex  integral,  hi  practice,  integral  methods  use 

3, 

some  convenient  approximation  based  on  Equation  7. 


A  detailed  critique  of  all  the  previously  reported  kinetic  treaments  is 
beyond  the  scope  of  this  report.  Several  reviews  have  been  published  (References 
2,  8,  and  9)  and  the  recent  paper  by  Flynn  and  Wall  (Reference  6)  offers  an 
excellent  critical  survey  of  many  of  the  methods.  It  is  fitting,  however,  to 
examine  some  of  the  approaches  and  evaluate  their  usefulness  for  obtaining 
meaningful  information  about  the  degradation.  Some  representative  approaches 
are  listed  in  Table  L  It  is  found  that  almost  all  the  methods  are  subject  to  one 
or  more  of  the  following  criticisms; 

a.  The  equation  fits  a  limited  portion  of  the  curve  so  that  parameters  for 
the  whole  degradation  range  cannot  be  obtained. 
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b.  Degradation  kinetics  are  assumed  to  be  of  the  “order”  type,  and  per¬ 
haps  even  first  order  kinetics  are  specified.  Since,  as  has  been  previously 
discussed,  this  assumption  is  not  necessarily  valid,  it  is  impossible  to  conclude 
whether  the  resulting  parameters  are  meaningful.  Often  a  random  type  deg¬ 
radation  and  probably  other  types  will  give  a  deceptively  good  fit  to  the  first 
order  approximation  over  most  of  the  conversion  range.  Errors  in  assumed 
order,  however,  give  vast  errors  in  the  other  kinetic  parameters;  thus,  in 
Reference  6,  a  random  curve  corresponding  to  an  activation  energy  of  60 
kcal-mole”^  was  analyzed  by  a  first  order  treatment  to  have  an  activation  energy 
ranging  from  66  to  175  kcal-mole“^  depending  on  conversion. 

c.  The  time  saving  advantage  of  single  point  and  single  curve  methods  is 
obvious  but  it  cannot  be  overemphasized  that  such  analyses,  when  applied  to  a 
new  system  of  unknown  kinetics,  may  be  absolutely  worthless  with  regard  to 
significant  kinetic  parameters. 

d.  Care  must  be  taken  when  deriving  approximate  solutions  since  a  sequence 
of  steps  involving  successive  approximations  may  give  a  final  working  equation 
which  is  not  particularly  applicable.  Assumptions  and  approximations  should 
be  kept  to  a  minimum. 

On  this  basis,  methods  involving  a  single  thermogram  can  be  discarded  and 
the  potential  of  methods  using  several  thermograms  examined.  The  method  of 
Reich  using  two  heating  rates  (Table  I)  is  unsatisfactory  since  no  information 
about  the  rate  law  can  be  obtained.  Chatter jee’s  method  involving  different 
sample  wei^ts  must  also  be  discarded  since,  in  addition  to  the  assumption  of 
order,  the  reasoning  is  erroneous  (Reference  37).  When  wei^t  terms  are 
e;q)ressed  in  the  proper  fractional  units,  his  treatment  collapses  completely,  hi 
the  equation  for  order  (Table  I),  it  can  be  visualized  that  “n”  would  always  be 
unity,  irrespective  of  the  mechanism  of  degradation,  excepting  the  possible 
influence  of  eiqierimental  irreproducibility  and  diffusion  controlled  wei^t  loss. 

Therefore,  the  choice  of  methods  is  reduced  to  two,  one  integral  and  the 
other  differential;  both  require  a  series  of  thermograms  to  be  obtained  using 
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different  heating  rates  and  appear  to  be  similar  in  potential  and  general  appli¬ 
cability.  This  also  represents  the  conclusion  of  the  Flynn  and  Wall  survey.  In 
this  report,  selectivity  was  carried  one  stage  further  according  to  the  following 
arguments; 


a.  In  the  past,  integral  methods  have  had  the  advantage  that  tedious  manual 
derivations  of  rates  of  weight  loss  could  be  avoided.  Since  many  thermo¬ 
balances  now  incorporate  automatic  differentiation  which  allows  simultaneous 
recording  of  TGA  and  differential  thermogravimetry  (DTG)  cxirves  and  since, 
lacking  this  equipment,  rates  can  be  computed  from  weight-loss  data  (see 
Section  III),  differential  methods  need  no  longer  be  avoided. 

b.  The  integral  treatment  involves  more  approximations  than  the  differential 
method.  Thus,  the  complex  integral  is  normally  approximated  to  a  convenient 
working  equation. 


c.  The  integral  treatments  assume  a  linear  heating  rate  throughout  the 
degradation.  This  assumption  is  not  strictly  valid  since,  for  example,  volatil¬ 
ization  of  gaseous  products  produces  a  cooling  effect  on  the  polymer  and  a 
lowering  of  the  heating  rate. 

d.  In  order  to  determine  the  rate  law  dependence,  F(W),  in  the  integral 
method,  experimental  curves  must  be  fitted  to  theoretical  curves  for  known 
kinetic  schemes.  Experimental  errors  and  kinetic  irregularities  (the  latter  are 
often  observed  during  the  initial  and  final  stages  of  weight  loss)  would  create 
difficulties  in  curve  fitting.  A  more  serious  disadvantage  is  that  the  true 
kinetic  scheme  might  be  different  from  the  available  theoretical  ones,  e.g.,  the 
isothermal  rate  curve  maximum  at  40%  conversion  for  polystyrene  (Reference  22). 
hi  the  study  of  the  complex  condensation  polymers  currently  being  used  for  higti 
temperature  applications,  new  types  of  rate  dependencies  may  well  appear  and 
would  necessitate  laborious  trial  and  error  curve  fitting  in  the  integral  method. 
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3.  “FRIEDMAN”  METHOD  OF  ANALYSIS 

The  foregoing  arguments  resxilted  in  the  selection  of  the  method  first  used 
by  Friedman  (Reference  15).  The  basic  equation  used  is  shown  in  Table  I  and  is 
the  logarithmic  e3q)ression  of  Equation  3 


log  [ 


dW 


(W^-Wf)  dt 


1  =  log  A  +  log  F(W)  - 


2.303RT, 


(8) 


B 


Equation  8  denotes  that  both  the  rate  and  temperature  corresponding  to  any 

specific  value  of  F(W)  are  dependent  on  the  heating  rate  employed.  If  it  is 

assumed,  for  the  moment,  that  the  chemistry  of  the  process  is  independent  of 

W  -W 

temperature,  then  F(W)  will  be  constant  at  any  weight  loss,  ^ 

of  heating  rate.  Thus,  by  measvtring  or  calculating  the  rates  and  temperatures 
corresponding  to  the  same  fractional  wei^t  loss  in  a  series  of  thermograms 
obtained  at  different  heating  rates  and  by  plotting  log(rate)  against  reciprocal 


temperature,  the  activation  energy  can  be  calculated  for  any  conversion.  An 
examination  of  the  changes  in  calculated  activation  energies  over  the  whole 
conversion  range  reveals  whether  the  kinetics  vary  with  conversion  although  it 
may  be  difficult  to  resolve  experimental  errors  and  real  changes  in  activation 
energy.  If  the  activation  energy  does  not  vary  much  over  the  entire  range,  an 
average  value,  E  ,  is  calculated.  Hence,  by  rearranging  Equation  8 


logAF(W)  =  log  A  +  log  F(W)  =log(rate)g  + 


2.303RT, 


B 


(9) 


Values  of  log  A  F(W)  are  calculated  at  each  conversion  for  the  varioxis  heating 
rates.  Theoretically,  there  should  be  no  heating  rate  variation,  A  and  F(W)  being 
considered  to  be  independent  of  temperature  in  the  simple  case,  but  ejqjerimental 
errors  usually  result  in  a  small  spread  of  log  A  F(W)  values  so  that  an  average 
value  is  again  calculated. 


against  log  (or  log  ttt — ^  for  degradations  producing  a  residue)  reveals 

o  o 


When  this  process  has  been  repeated  for  all  conversions,  a  plot  of  log  A  F(W) 

W  W-Wf 

-  (or  log  ^ 
o  'o 

the  rate  law.  In  the  programmed  TGA  method,  this  plot  removes  the  effect  of 
heating  rate  and  changing  temperature;  for  diagnostic  purposes,  it  is  analogous 
to  the  rate- conversion  curve  in  isothermal  treatments,  hi  many  cases,  therefore, 
the  appropriate  weight  functionality  will  be  elucidated  (Table  II). 
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TABLE  II 


F(W) 


EFFECT  OF  F(W)  ON  LOG  A  F{W)  PLOT 

LOG  A  F(W)  PLOT 


W  W  -W 

’  W  -W. 
o  of 


straight  line,  slope  =  1 


/  W  -W 

\  W-Wf  I 
o  I 


straight  line,  slope  =  n 


“random”  function 


curve,  maximum  at  log(0,74) 
(26%  conversion) 


Other  types  of  kinetics,  such  as  those  giving  a  rate  maximum  at  conversions 
other  than  26%,  would  be  just  as  easily  resolved  by  thiS' plot.  Kinetic  irregu¬ 
larities  are  revealed;  for  example,  early  wei^t  loss  caused  by  lower  activation 
energy  processes  normally  shows  as  a  steeply  falling  portion  in  the  initial 
stages  of  the  log  A  F(W)  cxirve.  In  fact,  any  true  change  in  the  kinetics  will 
produce  a  change  in  slope  of  the  plot  based  on  the  assumption  of  a  single 
activation  energy  process.  The  extreme  sensitivity  of  log  A  F(W)  to  changes 
in  E  explains  why  an  average  E  is  used  instead  of  individual  values  at  each 
conversion.  If  the  latter  were  used,  the  experimental  fluctuations  of  E^  would 
outweigh  the  effect  of  rate  on  log  A  F(W);  the  resulting  plot  would  be  very 
scattered  and  no  information  about  the  weight  functionality  woxild  be  obtained. 
Other  potential  sources  of  error  in  determining  the  form  of  F(W)  are  wrongly 
assigned  conversion  ranges  and  use  of  wrong  conversion  units.  These  may 
seem  obvious  but  the  former  can  be  easily  done  in  some  complex  thermograms 
and  the  latter  follows  from  perusal  of  the  variety  of  expressions  used  in  the 

past.  For  example,  in  the  original  description  of  this  method,  Friedman  used  a 

W-W. 

different  e3q)ression  of  concentration,  (  - =-),from  that  proposed  in  this 

^o 

report.  The  result  was  a  very  high  apparent  order  of  reaction. 

The  pre-e3q)onential  factor  can  be  easily  calculated  by  subtracting  the 
F(W)  value  from  the  log  A  F{W)  term.  The  determination  of  A  could  be  ac¬ 
complished  graphically  by  replotting  log  A  F(W)  against  log  F(W)  and  ex¬ 
trapolating  the  resulting  straight  line  to  log  F(W)  =  O.  The  ability  of  the  method 
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to  completely  analyze  more  complex  degradation  systems  with  some  typical 
complicating  features  is  discussed  in  the  following  paragraphs. 


a.  Activation  energy  changes  with  conversion 

Certain  degradation  mechanisms  involve  real  activation  energy  changes 
with  conversion  (References  3  and  5).  The  reason  may  be  a  dependence  of 
on  the  molecular  weight  which  itself  varies  with  conversion,  a  change  in  mech¬ 
anism  as  in  poly  (methyl  methacrylate),  or  a  change  in  structure  of  the  polymer. 
Suppose  a  smooth  increase  in  is  observed  in  the  E^-conversion  plot  obtained 
by  the  preceding  analysis.  The  assumption  of  an  average  E^^  to  calculate  log 
A  F(W)  would  invalidate  the  rate  law  determination.  Features  such  as  rate 
maxima  would  still  be  discerned  but  the  overall  curve  would  be  skewed  relative 
to  the  “theoretical”  curve.  Then,  it  is  conceivable  that  a  smooth  E^  profile 
might  be  tised  instead  of  an  average  E^  to  give  more  meaningful  information 
about  F(W). 

b.  Random  Degradations 

The  complex  functionality  of  weight  for  this  mechanism  (Equation  4)  casts 
doubt  on  the  validity  of  the  assvimption  in  the  kinetic  analysis  that  F(W)  is 
independent  of  temperatxire.  The  source  of  this  contention  is  that  a,  the  degree 
of  degradation  inherent  in  the  function,  contains  a  temperature  term.  It  can  be 
shown,  however,  that  a  is  independent  of  heating  rate  and  depends  only  on  the 
conversion.  Under  dynamic  conditions 

t 

-  J  kdt 

a  =  I  -  e  o 

T  E 

.AT 

a  =  I  -  e  ®  (lO) 

’o 

That  a  is  independent  of  B  follows  from  integral  treatments,  particularly  that 
in  Reference  21.  The  method  is  still  valid  therefore  for  the  random  case.  The 
derivation  of  a  pre-ejq)onential  factor  in  the  random  case  may  be  difficult  since 
the  various  parameters  needed  to  calculate  F  (W)  may  not  be  available.  Equation  4 
is  an  approximate  solution  dependent  on  certain  boundary  conditions  and 
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may  not  be  applicable  to  some  real  cases  althou^  the  overall  random  curve 
is  still  observed.  The  normal  method  of  obtaining  A  from 


Vdkt  / 


=  L/e  (  Reference  5) 


max 


and  in  this  treatment, 


[  log  A  F(W)]  =  logA  +  log 

•'max  '  ' 


(ID 


(12) 


must  be  used  cautiously. 


c.  Complex  Mechanisms 

The  treatment  of  two  of  the  more  straightforward  complex  cases  was  dis- 
cxissed  in  Reference  6.  The  first  case  involved  competitive  reactions  in  which 
the  rate  curve  and  thermogram  appeared  similar  to  that  for  a  simple  reaction 
except  for  irregular  trends  in  the  maximum  rate.  The  ability  of  the  method  to 
resolve  the  two  reactions  depends  on  how  different  the  individual  parameters 
are.  If  they  have  similar  orders  and  activation  energies,  it  is  do\abtfxil  whether 
any  resolution  could  be  achieved.  The  second  case  consisted  of  two  independent 
reactions,  each  of  which  could  be  observed  in  thermograms  obtained  at  low 
heating  rates.  The  corresponding  activation  energies  were  obtained  at  low  and 
high  conversions.  Only  the  two  methods  involving  several  heating  rates  showed 
any  success  in  resolving  these  cases.  Inreal polymer  degradations,  the  following 
complex  cases  have  been  observed  in  previous  and  current  work: 


a.  The  thermogram  consists  of  several  consecutive  steps  with  distinct 
plateaus  between  the  decompositions.  This  case  is  easily  dealt  with  by  treating 
each  step  individually  as  a  simple  case. 

b.  The  thermogram  exhibits  overlapping  reactions  and  the  DTG  curve  has 
several  maxima.  This  could  be  visualized  as  Case  a.  in  which  the  second  step 
commences  before  the  first  reaction  is  complete.  This  is  not  strictly  analogoxis 
to  the  “independent  reaction”  discussed  in  “Complex  Mechanisms”  since  the 
amount  of  each  reaction  may  depend  on  the  heating  rate.  In  such  a  case,  a  proper 
analysis  may  be  thwarted  although  relevant  information  would  still  be  obtained 
for  the  low  temperature  reaction  by  studying  the  initial  portion  of  the  weight  loss. 
An  example  of  this  type  is  presented  in  Section  VL 
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c.  A  complex  curve  somewhat  similar  to  b.  has  been  observed  for  some 
aromatic  polyesters  (Reference  23).  For  these  polymers,  the  major  weight-loss 
reaction  changes  smoothly  into  a  slow  char-forming  reaction  which  gives  the 
rate  curve  a  long,  high  temperature  tail.  Once  again,  it  is  predicted  that  only 
limited  information  will  be  obtained  and  further  discussion  must  await  detailed 
examination  of  actual  examples. 

4.  SIGNIFICANCE  OF  KINETIC  PARAMETERS 

The  influence  of  experimental  variables  in  programmed  thermogravimetry 
and  their  effect  on  the  resulting  kinetic  parameters  has  been  the  subject  of 
several  reviews  (References  2,  8,  19,  24,  and  25)  and  has  convinced  some 
authors  that  the  parameters  are  purely  empirical.  However,  employing  careful 
techniques  and  strict  standardization,  many  of  the  sources  of  error  such  as 
weiring  errors,  diffusion  effects,  and  differences  dependent  on  the  physical 
form  of  the  sample  can  be  minimized  or  eliminated.  The  sources  of  error,  which 
could  be  considered  appropriate  to  programmed  methods  and  not  to  isothermal 
techniques,  are  heating  rate  and  temperature  errors.  The  latter  can  be  removed 
by  good  experimental  procedure  and,  as  has  been  shown  previously,  the  former 
effect  is  removed  in  the  ultimate  analysis  by  the  method  chosen.  It  is  concluded, 
therefore,  that  the  programmed  method  used  should  give  information  equivalent 
to  that  obtained  isothermally  and  it  is  contended  that  this  information,  especially 
the  activation  energy,  can  be  meaningful  with  respect  to  stability  and  mechanism. 
Thus,  as  in  Reference  26,  overall  activation  energies  have  been  related  suc¬ 
cessfully  to  the  energies  of  individual  steps  comprising  the  reaction.  A  better 
proof  must  be  to  compare  the  experimental  activation  energy  to  that  observed  in 
conventional  kinetic  analysis.  This  should  be  feasible  when  a  polymer  and  its 
model  degrade  by  exactly  similar  mechanisms.  The  field  of  condensation 
polymers  may  contain  examples  satisfying  this  condition. 
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SECTION  m 

EXPERIMENTAL  PROCEDURE 

1.  POLYMER  SAMPLES 

Descriptions  of  the  methods  employed  for  the  preparation  and  purification 
of  samples  are  given  in  the  reports  which  describe  in  detail  the  results  for 
those  polymers. 

2.  APPARATUS 

The  thermobalances  used  were  the  Ainsworth  Models  AV  and  RV  which  gave 
fiill  scale  recorder  deflections  of  100  mg  and  10  mg,  respectively,  hi  most 
cases,  several  deflections  of  the  recorder  pen  were  necessary  to  follow  the 
complete  weight  loss  of  samples. 

A  sectional  diagram  of  the  thermobalance,  degradation  tube,  and  furnace 
is  given  in  Reference  27. 

The  temperature  programmer  used  was  the  West  Gardsman  Model  JGB 
Program  Controller  which  operated  a  proportioning  power  simply.  The  temper¬ 
ature  set  point  was  driven  by  a  cam  cut  to  give  close  to  a  linear  increase  in 
temperature  with  time.  Variation  in  program  rate  was  effected  by  changing 
gears  in  the  motor  to  cam  gear  train.  Some  of  the  approximate  program  rates 
selectable  were  75,  90,  150,  280,  and  450°  per  hour  but  the  actual  program 
rates  were  calculated  for  each  run  from  the  temperature-time  data. 

3.  TEMPERATURE  MEASUREMENT 

The  measurement  of  the  actual  temperature  of  a  material  undergoing  weight 
loss  presents  many  difficulties.  In  theory,  an  ideal  method  for  measuring  the 
temperature  would  be  to  surround  a  thermocoiq^le  bead  completely  with  the 
sample  and  measure  the  thermocouple  millivolt  output.  In  practice,  difficulties 
arise  especially  if  complete  loss  of  material  takes  place  during  degradation. 
In  this  case,  the  thermocouple  bead  becomes  more  and  more  exposed  to  the  heat 
source  as  weight  loss  occurs.  Further,  it  may  not  be  safe  to  assume  a  correct 
temperature  will  be  determined  even  when  the  sensor  is  completely  surrounded. 
This  will  be  dependent  upon  the  spectral  characteristics  of  the  sample  and  the 
heating  method  employed  (UV  or  IR). 


16 


AFML-TR-68-181 
Part  II 


Another  difficulty  is  involved  in  the  simultaneous  determination  of  weight 
and  temperature.  Torsion  of  wires  from  the  thermocouple  to  a  stationary  support 
will  alter  the  mass  reading  or  may  cause  noise  in  the  weight  record.  Methods 
differing  in  complexity  have  been  devised  to  overcome  these  defects  but  none  is 
entirely  satisfactory  (References  19  and  28). 

For  this  work,  it  was  decided  to  measure  temperatures  by  placing  a  thermo¬ 
couple  in  a  thermowell  as  close  to  the  sample  as  possible.  For  several  of  the 
pol3miers,  a  series  of  calibration  degradations  was  run.  Sample  temperature, 
measured  by  a  thermocouple  in  direct  contact  with  the  polymer,  was  recorded 
and  corrections  to  the  thermowell  temperatures  were  obtained. 

In  most  cases,  the  temperature  correction,  AT,  was  of  the  form 

AT  =  C  +  A/3 


A  and  C  are  constant 
/3  is  the  heating  rate 

The  temperature  corrections  were  usually  in  the  range  of  5  to  15°C.  Similar 
lags  have  previously  been  reported  (Reference  24). 


In  the  range  of  temperatures  over  which  weightless  occurs,  severe  temper¬ 
ature  lags  may  occur  (Reference  25)  since  large  quantities  of  heat  are  called 
for  during  an  endothermic  process. 


It  is  realized  that  temperatures  measured  in  this  work  are  only  approximate 
but  it  is  felt  that  by  standardizing  conditions  (sample  size,  crucible  and  furnace 
geometry,  etc.)  and  making  the  corrections  described,  errors  are  minimized. 
It  is  hoped  eventually  to  be  able  to  recalculate  this  data  making  corrections  for 
the  lag  during  the  endothermic  weight  loss.  This  must  await  the  relevant  ex¬ 
perimental  data. 
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4,  PROCEDURE 

The  sample  (usually  100  mg)  was  weighed  into  a  small  quartz  crucible  which 
was  then  suspended  in  a  quartz  degradation  tube  by  a  fine  nichrome  chain 
connected  to  the  balance  beam.  A  counterbalance  was  applied  to  the  opposite 
side  of  the  beam  making  sure  that  weigjits  at  least  equal  to  the  expected  weight 
change  were  suspended  on  the  beam.  After  the  apparatus  had  been  pumped  down 
to  a  pressure  below  0.1  micron  of  mercury,  the  furnace  which  surrounds  the 
degradation  tube  and  the  programmer  were  switched  on.  After  the  wei^t  change 
had  occurred,  programming  was  continued  until  a  good  final  weight  base  line 
was  recorded.  From  the  thermogram  which  recorded  both  weight  and  temperature 
as  a  function  of  time,  the  rates  of  weight  loss  as  a  function  of  the  instantaneous 
percent  wei^t  loss  were  computed  (Section  IV).  hi  any  cases  in  which  a  steady 
final  weight  line  could  not  be  obtained  (c.g.i  for  some  aromatic  polyamides 
and  polyesters),  rates  were  based  on  the  initial  sample  weight  instead  of  on  the 
total  weight  loss. 

Occasionally  slight  initial  weight  losses  were  noticed  due  to  removal  of 
solvents  or  water  from  the  polymers.  In  these  cases,  data  were  taken  from  the 
thermogram  after  the  weight  line  was  again  level. 

In  an  attempt  to  reduce  procedural  errors,  the  ejq)erimental  procedure  and 
apparatus  were  standardized  as  far  as  possible  (e.g.,  furnace,  method  of  tem¬ 
perature  measurement,  sample  size,  crucible  geometry). 

Extremely  useful  information  on  experimental  methods  and  apparatus  iised 
in  thermogravimetry  is  given  in  Reference  25. 
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SECTION  IV 

COMPUTER  HANDLING  OF  THERMOGRAVIMETRIC 

ANALYSIS  DATA 


This  section  is  devoted  to  the  treatment  of  TGA  data  using  Friedman’s 
method  which  is  discussed  at  length  in  Section  II.  There  is  some  similarity  to 
the  treatment  of  isothermal  thermogravimetry  data  (Reference  35).  Despite 
the  possible  repetition,  this  section  describes  all  aspects  of  the  handling  of 
programmed  thermogravimetry  data.  The  first  of  this  section  is  concerned  with 
the  determination  of  the  rates  of  wei^t  loss  and  the  second  portion  describes 
the  evaluation  of  the  parameters  involved  in  Equation  9. 

1.  COMPUTATION  OF  RATE  OF  WEIGHT  LOSS 

hi  Section  m  the  experimental  procedure  is  described.  Figure  2  shows  a 
typical  recorder  trace  from  the  thermobalance.  The  two  curves  represent 
temperature  (measured  by  a  Chromel/Alumel  thermocouple  located  in  a 
thermowell  close  to  the  sample  container)  and  the  sample  weight  (measured 
electronically  by  determining  changes  in  the  resonant  frequency  of  a  transducer 
caused  by  deflection  of  the  balance  beam).  The  pen  excursions  are  linearly 
dependent  i^ion  temperature  and  weight,  full  scale  deflections  corresponding 
to  500  or  lOOO’C  and  10  mg  or  100  mg  weight  change. 

Obviously  the  two  pens  cannot  travel  on  the  same  line  perpendicular  to  the 
time  axis.  A  small  correction  has  to  be  made  to  data  read  from  the  same  line 
to  ensure  that  pairs  of  weights  and  temperature  data  represent  conditions  at 
the  same  time. 


Methods  are  available  for  automatically  converting  signals  from  measuring 
equipment  to  digital  form  for  computer  processing,  but  such  methods  were  not 
on  hand  for  this  work.  Thus  it  was  necessary  to  obtain  the  recorder  traces  and 
to  take  data  from  the  two  curves  either  using  a  mechanical  graph  reader  or 
manually,  hi  the  latter  method,  the  chart  was  taped  to  a  board  and  scales 
graduated  in  suitable  increments  (20th  or  32nd  of  an  inch)  taped  to  both  sides  of 
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the  chart  parallel  to  the  time  axis.  A  ruler  was  then  moved  along  the  time  axis 
and  corresponding  wei^t  and  temperature  data  recorded  at  200  to  300  positions, 
hi  order  to  simplify  data  gathering,  elapsed  inches  of  chart  were  recorded  and 
time  regenerated  during  the  machine  calculations  using  the  chart  speed  (in  inches 
per  minute).  The  zero  time  data  is  taken  from  the  curve  at  some  suitable 
position  prior  to  significant  weight  loss. 

The  first  card  identifies  the  deck  of  cards  that  follows.  The  information 
and  the  columns  in  which  it  is  punched  is  shown  in  Table  IIL 

TABLE  in 

IDENTIFICATION  CARD 


COLUMNS 

CONTENTS 

1-8 

Any  identification,  name,  etc. 

9-12 

Run  number  (e.g.,  04P) 

13  -  14 

Blanks 

15  -  22 

Date 

23  -  24 

Blanks 

25  -  45 

Alphabetic  information,  polymer  name,  etc 

46  -  52 

Chart  speed  in  inches  per  minute 

53  -  54 

Blanks 

55  -  59 

Pen  offset  in  inches 

60  -  62 

hiches  between  first  and  last  data 

63  -  70 

Blanks 

71  -  73 

Temperature  at  initial  time  reading 

74  -  76 

Blanks 

77  -  79 

Temperature  at  final  time  reading 
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Following  the  identification  card  is  the  data  deck.  The  data  cards  are 
punched  as  shown  in  Table  IV.  A  typical  input  card  deck  is  reproduced  in 
Appendix  I. 


TABLE  IV 

CONTENTS  OF  DATA  CARDS 


COLUMNS 

CONTENTS 

1-13 

Identification  information,  name,  etc. 

14 

An  index,  LBJ,  to  signal  the  end  of  a  deck 

15  -  17 

Blanks 

18  -  22 

Time  in  inches 

23  -  28 

Weight 

29  -  32 

Temperature 

33 

Blank 

34  -  38 

Time  in  inches 

39  -  44 

Wei^t 

45  -  48 

Temperature 

49 

Blank 

50  -  54 

Time  in  inches 

55  -  60 

Wei^t 

61  -  64 

Temperature 

65 

Blank 

66  -  70 

Time  in  inches 

71  -  76 

Weight 

77  -  80 

Temperature 

Using  the  input  data,  the  computer  assembles  a  matrix  of  time,  wei^t, 
and  temperature  data  which  is  then  scanned  to  find  the  next  weight  after  a  1% 
weight  loss.  A  niunber  of  pairs  of  weight  and  time  data  either  side  of  this  point 
is  then  fitted  to  a  quadratic  (using  PLSQ  least  squares  curve  fit  subroutine).  The 
quadratic  is  then  solved  for  the  time  taken  for  an  exact  1%  weight  loss  to  occur. 
These  calculations  are  repeated  for  all  integral  percent  weight  losses  xq)  to  99. 
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Time  and  temperatxire  data  are  fitted  to  a  polynomial  using  the  PLSQ 
subroutine.  A  single  high  order  polynomial  is  used  to  fit  all  the  temperature¬ 
time  data  whereas  a  quadratic  is  used  to  fit  short  sections  of  the  weigfit-time 
data.  These  two  different  methods  were  chosen  because  a  quadratic  may  easily 
be  solved  whereas  hi^er  orders  may  not.  Only  substitution  into  the  time- 
temperature  polynomial  is  needed  here  so  high  orders  may  be  used  to  get  a 
better  fit  of  the  data.  The  weight-time  curve  fit,  however,  had  to  be  solved  to 
interpolate  time  values;  as  a  quadratic  fit  was  dictated,  the  best  fit  was  achieved 
by  using  a  relatively  small  number  of  curve  fit  data. 

After  the  input  data  has  been  compared  with  results  calculated  from  the 
fitted  curves,  the  spurious  input  is  replaced  by  fitted  values  and  the  curve  fits 
recalcTilated.  The  rate  of  weight  loss  is  then  calculated  for  each  percent  wei^t 
loss  and  the  average  heating  rate  computed.  A  schematic  representation  of  this 
rate  program  is  given  in  Figure  3,  and  the  complete  program  is  reproduced  in 
Appendix  IL 

Finally  the  results  are  printed  out  in  tabular  form  together  with  graphs 
representing  the  variations  of  rate  of  wei^t  loss  with  percent  weight  loss, 
weight  with  time,  and  temperatxire  derivative  with  time.  Results  are  also  punched 
on  to  computer  cards  for  further  processing.  A  typical  punched  card  output 
deck  is  shown  in  Appendix  III  and  examples  of  pointout  data  is  given  in 
Appendix  V. 

2.  COMPUTATION  OF  ARRHENIUS  PARAMETERS 

Friedman’s  method  for  calculating  activation  energy,  etc.,  which  is  de¬ 
scribed  fully  in  Section  II,  requires  the  comparison  of  wei^t-loss  rates  ob¬ 
tained  from  a  number  of  thermogravimetric  analyses  carried  out  at  differing 
rates  of  temperature  increase.  Comparison  is  made  between  rates  of  wei^t 
loss  at  identical  extents  of  weight  loss,  and  activation  energy  may  then  be 
calculated  from  the  slope  of  the  curve  of  log  (rate  of  weight  loss)  against 
reciprocal  of  the  absolute  temperature.  The  lower  the  heating  rate,  the  lower 
the  temperature  will  be  for  a  given  wei^t  loss. 
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Figure  3,  Schematic  Representation  of  Rate  of  Weight  Loss  Calculations 
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The  output  punched  cards  from  the  rate  program  contain  both  rate  of 
weight  loss  and  temperature  for  each  percent  wei^t  loss  so  the  program  merely 
selects  corresponding  data  from  each  experiment  and  carries  out  a  least  squares 
straight  line  fit  of  the  log  (rate  of  weight  loss)  and  l/T  data.  Both  the  slope  and 
intercept  of  the  best  line  are  computed.  Log  A  F(W)  values  are  then  calculated 
using  an  average  value  for  the  activation  energy.  The  range  of  weight  loss  over 
which  this  average  is  computed  is  usually  chosen  to  omit  very  low  and  high 
conversions.  A  discussion  of  the  effects  of  changes  in  activation  energy  with 
conversion  is  in  Section  IL 

Figure  4  shows  a  schematic  representation  of  Arrhenius  parameter  cal¬ 
culations,  and  in  Appendix  IV  the  complete  Arrhenius  program  is  reproduced. 
Typical  resiilts  from  all  these  calculations  are  presented  and  discussed  in 
Section  V. 
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Read  In  Identification  Card 


Read  in  Data  Card 


Check  for  Consectutive  Order  of  Input 


Decide  If  Lott  Card  Has  Been  Reached 


Check  for  any  Zero  Rates.  If 
Present,  Set  Up  Dummy  Points 
for  Graph  Plot  (GP) 


Compute  Log  (Rate  of  Weight  Loss) 
and  Reciprocal  of  Absolute 
Temperature 


Least  Squares  Determination  of 
Slope  and  Intercept  of  Plot  of 
Log  (Rote  of  Weight  Loss) 
Against  Reciprocal  of  Temperature 


Calculate  Average  Activation 
Energy  and  Pre-exponential  Factor 


Calculate  A  F(ai)Ueing  Average 
Activation  Energy 


Write  Out  Results,  Percent  Weight 
Loss,  Activation  Energy,  Pre-exponential 
Factor,  Average  A  F(ai)and 
Individual  A  F(a;)Values 


Plot  Graphs  of  Log  (Rate  of  Weight  Loss) 
Against  Reciprocal  of  Temperature, 
Activation  Energy  and  Pre-exponentiol 
Factor  Against  Percent  Weight  Loss, 
Average  A  F(cu) Against  Log 
(Percent  Weight  of  Residue) 


Figure  4.  Schematic  Representation  of  Arrhenius  Parameter  Calculations 
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SECTION  V 

WEIGHT  LOSS  OF  POLY(TETRAFLUOROETHYLENE) 

To  validate  the  procedure  for  the  determination  of  activation  energy  and 
order  of  reaction  which  is  described  in  Section  II,  results  obtained  by  this 
method  were  compared  with  those  previously  published  for  the  degradation  of 
poly(tetrafluoroethylene),  Teflon.  This  particular  compound  was  chosen  for  its 
relative  lack  of  complications  and  for  the  availability  of  published  information. 

Madorsky  and  co-workers  (Reference  29)  measured  the  weight  of  samples 
of  Teflon  maintained  at  fixed  temperatures.  Plots  of  the  rate  of  wei^t  loss 
against  the  percent  volatilization  were  linear  between  about  20  and  80%  vola¬ 
tilization  showing  the  degradation  to  be  a  first  order  process.  Their  kinetic 
data  could  be  summarized  by  the  following  equation 


k 


I 


4.7  X 


I0'®e 


80,500/RT  -1 

’  sec 


k^  =  first  order  rate  constant 
R  =  gas  constant 
T  =  absolute  temperature 

A  mechanism  for  degradation  involving  thermal,  weak  link,  or  end  initiation 
followed  by  unzipping  of  the  free  radicals  produced  was  postulated. 

Wall  and  Michaelsen  (Reference  30)  confirmed  these  observations  but  sug¬ 
gested  that  a  zero  order  dependence  of  the  rate  of  weight  loss  on  sample  weight 
was  observable  below  about  480‘’C.  They  presented  data  which  showed  that  at 
460‘’C  under  nitrogen  the  weight  loss  of  Teflon  is  a  linear  function  of  time  up  to 
about  40%  weight  loss. 


Anderson  (Reference  31)  analyzed  Teflon  thermogravimetry  data  by  the 
method  of  Freeman  and  Carroll  (Reference  11)  and  found  the  degradation,  in 
vacuum,  to  be  first  order  between  450  and  550°C  with  an  activation  energy  of 
75  ±4  kcal/mole. 
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Reich  and  co-workers  (Reference  13)  analyzed  Teflon  weight-loss  data 
using  the  method  of  Anderson  and  Freeman  (Reference  12)  previously  discussed. 
They  obtained  activation  energies  varying  between  69  and  74  kcal/mole,  the 
average  being  72  kcal/mole. 

Lee  and  co-workers  (Reference  32)  presented  information  obtained  using 
heating  rates  between  300  and  1200'’C  per  hour.  It  wiU  be  shown  here  that  such 
hi^  heating  rates  are  likely  to  result  in  large  uncertainties  in  the  temperature 
measurement  with  consequent  curvature  of  activation  energy  plots.  However 
they  quote  between  60  and  69  kcal/mole  and  orders  of  reaction  between  0.7 
and  0.85  depending  on  the  method  of  plotting  employed. 

Carroll  and  Manche  (Reference  33)  re-examined  Madorsky’s  data  and  deter¬ 
mined  the  activation  energy  as  a  function  of  the  conversion  and  showed  that 
between  10  and  80%  weight  loss  the  activation  energy  decreases  from  80  to  about 
46  kcal/mole.  The  decrease  of  E  with  increasing  conversion  was  apparent  for 

cL 

both  the  programmed  temperature  increase  and  for  the  isothermal  weight  loss 
of  Teflon  in  vacuum.  The  reaction  was  said  to  be  zero  order. 

Section  II  gives  a  detailed  discussion  of  the  various  techniques  which  have 
been  used  for  calculating  kinetic  parameters  from  weight-loss  data.  Our  con¬ 
clusions  are  that  each  method  has  serious  shortcomings.  We,  therefore, 
consider  that  the  data  of  Madorsky  (Reference  29)  which  was  obtained  iso- 
thermally  should  be  the  most  reliable  for  comparison  with  the  results  of  our 
investigation. 

As  a  check  on  the  present  differential  method  for  the  determination  of  E 

a 

and  order  of  reaction,  a  series  of  programmed  temperature  increase,  vacuum 
weight  loss,  experiments  on  100  mg  samples  of  Teflon  was  carried  out  using 
heating  rates  between  45  and  450°  per  hour.  The  polymer  used  was  Du  Pont 
Teflon  molding  powder  Composition  6  in  the  form  of  fine  granules. 

The  first  runs  which  were  carried  out  at  the  hi^er  heating  rates  gave  very 
high  rates  of  weight  loss  (~  10%  per  minute).  With  such  high  rates  of  reaction, 
considerable  temperature  lags  migiit  be  expected  so  measurements  were  also 
made  using  very  low  heating  rates  (below  150°  per  hour). 
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Figure  5  is  a  plot  of  log  (rate  of  weight  loss)  against  the  reciprocal  of  the 
absolute  temperature  for  50%  conversion.  Plots  for  other  conversions  were 
similar  to  this  one.  It  is  evident  that  a  straight  line  cannot  be  drawn  to  represent 
the  data  over  the  whole  temperature  range.  However,  the  runs  carried  out  at 
150“  per  hour  and  lower  heating  rates  do  show  a  linear  dependence  of  log  (rate 
of  weight  loss)  on  l/T.  The  slope  of  the  line  drawn  through  only  these  four 
points  gives  an  activation  energy  of  69.3  kcal/mole.  A  case  could  be  made  for 
considering  only  the  three  lowest  heating  rates.  The  derived  activation  energy 
would  then  be  increased. 

It  is  probable  that  the  curvature  of  the  Arrhenius  plot  when  the  higher 
heating  rates  are  \ased  is  due  to  the  large  thermal  lags  when  the  rates  of  wei^t 
loss  are  large.  In  the  hope  of  bringing  these  results  into  line,  an  attempt  will  be 
made  to  correct  for  these  lags  by  direct  sample  measurement. 

The  computer  printout  for  the  rates  of  wei^t  loss  for  one  of  the  Teflon 
ejqieriments  is  given  in  Appendix  V  and  Appendix  VI  shows  the  rates  of  wei^t 
loss  at  each  1%  conversion  for  the  four  lowest  heating  rate  runs.  Appendix  VII 
is  the  computer  printout  for  the  activation  energy  calculations  based  i?)on  all 
the  rates  quoted  in  Appendix  VL  Figure  5  contains  all  the  50%  data  from 
Appendix  VI  as  well  as  data  obtained  using  hi^er  heating  rates. 

1.  VARIATION  OF  ACTIVATION  ENERGY  WITH  CONVERSION 

Figure  6  shows  a  plot  of  activation  energy  as  a  function  of  wei^t  loss,  the 
data  being  taken  from  Appendix  VII.  Between  4  and  99%  wei^t  loss,  activation 
energy  varies  between  about  62  and  83  kcal/mole.  The  average  value  of 
69.34  kcal/mole  between  10  and  80%  weight  loss  has  been  used  to  calculate 
log  A  F(W)  values  used  in  one  of  the  curves  in  Figure  7.  Careful  inspection 
of  the  activation  energy  data  reveals  an  approximately  constant  value, 
average  =  63.98  kcal/mole,  between  10  and  50%  weight  loss.  At  greater  con¬ 
versions,  E  increases  slowly  to  a  maximxun  which  is  maintained  between  65 

Si 

and  80%  weight  loss. 

It  has  been  reported  that  Teflon  undergoes  a  change  in  physical  properties 
above  50%  weight  loss.  In  Reference  30  it  is  claimed  that  the  polymer  melts 
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Figure  5.  Arrhenius  Plot  for  Degradation  of  Teflon  at  50%  Weight  Loss 
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Figure  6.  Variation  of  Activation  Energy  with  Conversion  for  the  Weight  Loss 
of  Teflon 


Log  AF(W) 


ArML-TR~68-181 
Part  II 


AFML-TR-68-181 
Part  n 

at  this  conversion.  There  may  be  some  correlation  between  this  phenomenon 
and  the  change  in  the  activation  energy  curve  at  the  same  conversion. 

2.  ORDER  OF  REACTION 

Figure  7  also  shows  the  log  A  F(W)  obtained  using  the  activation  energy 
average  for  the  10  to  50%  weight  loss  (i.e.,  63.98  kcal/mole).  The  straight  lines 
drawn  which  represent  the  20  to  80%  weight  loss  data  fairly  well  have  slopes  of 
0.84  and  0.89,  and  only  sligjit  maxima  in  the  log  A  F(W)  curves  are  apparent  at 
very  low  conversions.  Thus  the  weigjit-loss  process  obeys  approximately  first 
order  kinetics.  Other  lines  having  differing  slopes  may  be  drawn  throu^  points 
representing  more  limited  weight  loss  ranges  but  it  would  be  unwise  to  assign 
real  significance  to  such  slopes. 

3.  PRE-EXPONENTIAL  FACTOR 

The  extrapolated  value  of  log  A  F(W)  at  zero  percent  weight  loss  is  a 
measure  of  the  pre-exponential  factor  (A).  Using  the  low  E  value,  log  A  =  15,7, 

Si 

with  the  higfi  E^,  log  A  =  17.1,  (Note:  Since  weight  losses  used  here  are  in 
percent,  log  100  has  been  subtracted  from  each  intercept.) 

4.  DISCUSSION 

The  data  presented  here  is  in  fair  agreement  with  some  of  the  published 
information.  The  10  to  80%  average  activation  energy  is  about  10  kcal/mole 
lower  than  Mador sky’s  value.  The  change  in  E  with  conversion  does  not  agree 
with  the  drop  from  80  to  46  kcal/mole  calculated  by  Carroll  and  Manche 
(Reference  33).  It  is,  however,  more  reasonable  to  expect  the  increase  in  E 

Si 

with  conversion,  as  we  find,  then  to  e3q)ect  a  large  drop  in  E  . 

If  the  degradation  involves  unzipping  with  a  long  kinetic  chain  length 
throughout  the  total  weight  loss,  a  first  order  rate  dependency  would  be  indi¬ 
cated,  and  no  changes  in  molecular  weight  of  the  residue  would  occur.  It  has  been 
shown  (Reference  34)  that,  at  500®C,  tetrafluoroethylene  is  the  major  volatile 
product  of  degradation  (95%)  but  small  amounts  of  CF.  and  CoF„  are  also 
produced.  At  higher  temperatures,  the  yield  of  tetrafluoroethylene  is  reduced, 
other  products  being  produced  by  termination  of  short  kinetic  chain  length 
unzipping  processes.  If  the  kinetic  chain  length  is  shorter  than  the  degree 
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of  polymerization,  a  change  in  molecular  weight  of  the  residue  would  take  place 
during  weight  loss  with  consequent  complication  of  the  degradation  mechanism. 

It  is  felt  that  the  method  described  here  for  the  determination  of  kinetic 
parameters  involved  in  thermogravimetry  gives  adequate  agreement  with 
literature  data  for  Teflon  to  justify  its  application  to  other  systems.  The 
kinetics  of  degradation  of  pol3rmers  which  obey  more  complicated  laws  are 
disciissed  in  Section  VI. 


SECTION  VI 

DEGRADATION  OF  OTHER  POLYMERS 

In  this  section,  representative  examples  of  results  obtained  during  thermal 
degradation  of  polycondensates  are  discussed.  These  examples  have  been  chosen 
to  show  several  types  of  log  A  F(W)  curves  which  are  derived  from  weight- loss 
data  using  the  computational  methods  described  in  Section  IV. 

1.  POLY  (1,4-  PHENYLENE  SEBACATE) 

Figure  8  shows  how  the  rate  of  weight  loss  for  this  polymer  varies  with  the 
extent  of  conversion.  That  the  mechanism  of  the  degradation  of  this  polymer  is 
more  complex  than  that  for  poly(tetrafluoroethylene)  is  shown  by  the  fact  that 
the  curve  of  rate  of  weight  loss  against  percent  weight  loss  exhibits  two  distinct 
maxima,  one  at  about  45%  and  the  second  at  about  90%  of  the  overall  wei^t  loss. 
Since  separation  of  the  maxima  is  apparent,  the  activation  energies  of  the  in¬ 
dividual  processes  must  differ  appreciably.  The  greater  the  energy  difference, 
the  better  will  be  the  resolution  of  the  rates  of  each  process. 
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Percent  Weight  Loss 

Figure  8,  Rate  of  Weight  Loss  of  Poly(l,4-phenylene  sebacate)  as  a  Function 
of  Weight  Loss 
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Figure  9  shows  that  the  activation  energy  does  indeed  vary  during  the 
degradation  of  poly(l,4-phenylene  sebacate).  hi  the  early  stages  of  degradation 
(Line  A),  is  about  30  kcal/mole  and  a  gradual  rise  takes  place  until  is 
close  to  60  kcal/naole  during  the  final  20%  weight  loss  (Line  B). 

The  curve  of  log  A  F(W)  against  log  (residue  weigfit-percent),  produced 
when  the  overall  average  E  of  36.1  kcal/mole  is  used,  is  shown  as  Curve  III 

d. 

in  Figure  10.  A  change  in  the  slope  of  this  curve  is  apparent  at  conversions 
exceeding  about  70%  showing  the  change  in  mechanism  brought  about  by  the 
commencement  of  the  second  reaction.  As  explained  in  Section  II,  an  erroneously 
high  slope  would  be  derived  from  this  curve  since  this  type  of  plot  should  be 
based  on  the  weight  loss  during  a  single  component  reaction,  here  the  weight-loss 
data  is  based  on  the  sum  of  the  two  component  reactions. 

In  order  to  separate  the  contributions  due  to  each  of  these  two  processes, 
it  is  necessary  to  go  back  to  the  original  curve  of  rate  of  weight  loss  as  a 
function  of  overall  percent  conversion  (Figure  8)  or  the  curve  of  rate  of  weight 
loss  against  temperature  (Figure  11).  By  careful  inspection  and  judicious  use  of 
curve  drawing  techniques,  it  is  possible  to  resolve  the  two  peaks,  from  either 
curve,  into  the  pure  components.  Overlap  of  the  two  reactions  occurs  between 
about  45  and  80%  of  the  total  weight  loss.  When  a  separation  has  been  made, 
calculations  of  log  A  F(W)  for  each  component  may  be  made  using  the  relevant 
conversions.  The  required  activation  energies  are  found  from  the  approximately 
linear  parts  of  the  E  against  weight-loss  curve  (Lines  A  and  B  in  Figure  9). 

Figure  10  shows  the  results  of  such  a  resolution  of  a  complex  weight-loss 
process.  Curve  I  is  the  log  A  F(W)  plot  for  Component  Reaction  L  This  process 
is  probably  random,  the  low  conversion  rise  in  this  ciirve  probably  being  due 
to  an  early  low  activation  energy  weight-loss  process.  Cxmve  II  represents 
Component  Reaction  IL  The  drop  in  the  curve  at  low  conversion  may  not  be 
significant  since  this  is  the  region  of  maximum  overlap  with  Reaction  L  The 
slope  of  the  curve  at  higher  conversions  shows  the  reaction  obeys  either  first 
order  or  random  kinetics. 
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Figure  10.  Log  A  F(W)  Curves  for  Overall  and  Component  Reactions  in  the 
Degradation  of  Poly(l,4-phenylene  sebacate) 
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Figure  11,  Separation  of  Rate  of  Weight  Loss/Temperature  Curves 
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Thus  it  can  be  seen  that  it  is  possible  to  separate  to  some  extent  the  compo¬ 
nents  of  a  complex  weight-loss  process.  However,  where  overlap  occurs  there 
is  some  doubt  about  the  rates  of  the  individual  wei^t  losses.  This  doubt  is 
reflected  in  the  inability  to  assign  real  significance  to  the  part  of  the  log  A  F(W) 
curve  which  includes  considerable  overlap  of  reactions. 

A  case  could  be  made  from  the  data  shown  in  Figure  8  that  another  compo¬ 
nent  may  be  present  in  the  range  60  to  70%  weight  loss  but  even  if  real  it  would 
be  virtually  impossible  to  achieve  its  resolution. 

A  similar  treatment  of  the  same  data  is  given  in  Reference  23. 

2.  POLY(HEXAMETHYLENE  SEBACAMIDE),  NYLON  6.10 

Unlike  many  aromatic  polyamides  (Reference  40),’  nylon  6.10  degrades 
completely  leaving  no  appreciable  residue.  The  maximum  rate  of  weight  loss 
occurs  at  about  60%  weight  loss. 

Figure  12  shows  the  variation  of  activation  energy  with  the  extent  of  weight 
loss  determined  from  programmed  thermogravimetry  data.  The  activation 
energy  rises  rapidly  during  the  first  10%  weight  loss  and  then  remains  in  the 
region  of  57  kcal/mole  for  the  remainder  of  the  weight-loss  process.  The  early 
rise  in  E  can  probably  be  attributed  to  the  removal  of  absorbed  water  from  the 

d 

pol3aner  or  distillation  of  low  molecular  wei^t  volatiles.  Using  an  activation 
energy  of  57.2  kcal/mole,  the  log  A  F(W)  curve  shown  in  Figure  13  has  been 
constructed.  The  20  to  90%  wei^t-loss  data  is  represented  by  a  good  straight 
line  having  a  slope  of  1.02  indicative  of  a  random  or  first  order  decomposition 
mechanism.  The  downward  curvature  of  the  line  at  low  conversions  tends  to 
indicate  a  random  process  is  operative  especially  as  the  maximum  occurs  close 
to  25%  weight  loss  (Reference  5).  Other  workers  have  concluded  that  the  same 
mechanism  describes  the  degradation  of  other  polyamides  but  a  possible  ionic 
hydrolysis  process  may  occur  simultaneously  (Reference  38). 

Several  other  aliphatic  polymers  have  been  studied  and  found  to  yield 
similar  log  A  F(W)  curves  and  mechanistic  interpretations.  The  degradation  of 
several  aliphatic  polyesters  and  polyamides  is  discussed  in  References  36  and  39. 
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SECTION  vn 
CONCLUSIONS 


The  method  of  Friedman  has  been  shown  to  be  applicable  to  a  wide  variety 
of  polymer  degradations.  It  can  detect  a  single  first  order  wei^t  loss,  one 
where  random  decomposition  is  suspected  (as  noted  by  a  maximum  in  the  rate 
law  curve)  and  can  provide  some  resolution  of  a  wei^t-loss  curve  showing 
more  than  one  rate  maximmn.  It  is  our  conclusion  that  only  throu^  the  use  of 
several  thermograms  at  different  heating  rates  can  this  information  be  obtained 
reliably. 
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z 

c 

PROGRAMTO”  DETERMINE  RATFS  OF  WEIGHT  LOSS  AT  ONE  PERCENT  WEIGHT  LOSS 
INTERVALS  FROM  THERMOGRAV IMETR IC  ANALYSIS  DATA. 

t 

c 

INPUT  TEMPERATURES  ARE  FITTED  TO  A  FIFTH  DEGREE  POLYNOMIAL  USING  A 

LEAST  SQUARES  SUBROUTINE  (PLSQ). 

c 

c 

WEIGHTS  CORRESPONDING  TO  SHORT  TEMPERATURE  RANGtS  ARE  FITTED  TO  A 

QUADRATIC  BY  PLSQ. 

L 

C 

INPUT  WEIGHTS  DIFFERING  FROM  FITTED  LINE  BY  MORE  THAN  OnE  PERCENT  OF  THE - 

TOTAL  WEIGHT  LOSS  ARE  REPLACED  BY  THE  CURVE  FIT  VALUE. 

L 

C 

OUTPUT  DATA  IS  PUNCHED  ON  TO  CAROS  FOR  FURTHER  PROCESSING  (TO  CALCULATE 
ACTIVATION  ENERGY  ETC). 

L 

C 

INPUT  DATA.  SPEEO=CHART  SPEED  IN  INCHES  PER  MINUTE,  0FFS5T=PEN  SEPARATION 

IN  INCHES.  ITI=INITIAL  TEMPERATURE  READING,  ITF=FINAL  TEMP.  READING. 

L 

C 

INDEX  LBJ  IS  SETE(Ju4L  TO  1  TO^STCTN AL  TAST  CARD  UF  A~RI]N  '  ^  ^ 

L 

DIMENSION  T(1509) , W( 1500 ) ,TEMP(  1500  ) ,WW( 105 ), TNW(  105 ), DWDT( 105 ) , 

.X(5C),Y(50),C(10),A(l),PLQT(50,105),Z(50),TOER{105),TPOLYa05), 

.RTEMP ( 105) 

1 

READ  (5,1000)  ID, DATE1,DATE2,C0M1, COM2, COM3, CDM4,  SPEED , OFF SE T ,  '  ' 

.XRANGE , ITI  ,  ITF 

1 

CHART  =  2.545454*SPEED 

WRITE  (6,3000) 

5 

WRITE  (6,3050)  I D , D ATEl , D ATE2 , COMl, C0M2,C0M3, COM4 

WRITE  (6,3060)  SPEED 

6 

7 

WRITE  (6,3065)  CHART  . .  ’  . 

WRITE  (6,3070)  OFFSET 

8 

9 

WRITE  (6,3080)  XRANGE 

WRITE  (6,3090)  ITI 

10 

11 

WRITE  (6,3100)  ITF 

N  =  0 

12 

L 

C 

START  LOOP  TO  READ  IN  DATA  CARDS 

C 

DO  30  1=1,2995,4 

N  =  N  +  4 

READ  (5,1010)  LBJ,T(I ) ,W( I ) ,TEMP( I ) ,T( I+1),W( I+l) ,TEMP{ I+l) , 

C 

.  I  (  1+2)  ,W(  1  +  2)  ,TcMP(  1+2  )  ,T(  1+3)  ,W(  I-*-3),TtMP(  1  +  3)  . . .  '  ’  . . 

C 

C 

T  =  TIME  (IN  INCHES),  W  =  WEIGHT,  TEMP  =  TEMPERATURE  . . 

17 

30 

IF  (LBJ.EQ.l)  GO  TO  40 

CONTINUE 

C 

'40 

NL  =  N-4  . . . . .  " 

L 

C 

START  LOOP  TO  CALCULATE  NUMBER  OF  DATA  READ  IN 

DO  5C  I=NL,N 

IF  (T(I  ).NE.O.O)  GO  TO  50 

,M  =  i-i  --  "  -  -  - 

GO  TC  51 

50 

51 

CONTINUE 

JJ  =  .030*FL0AT(N) 

C 

'LL  =  MAX0(JJ,10) 

L 

JJ  =  3  PERCENT  OF  NO.  OF  DATA  SETS  READ  IN  -  -  . 

50 
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plot _ ^  EFN _ SOURCE  STATEMENT  -  IFIM(S)  - _ 

LL  =  NO.  OF  ^_RyE  f  POINTS  (LATER  =  NN) _ 

WRITE  (6,3170)  LL  _  _ 

WRITE  {6,3010)  N 

NN  =  LL  _ 

K  =  POLYNOMAL  ORDERj,  NEED^p^OR  Plj^Q^  ^BjURTT  I N  LIST  =  0  FOR  NO  ERROR _ 

ANALYSrS  OF  PLSG 

D  =  total  W_EIGHT  LOSS _ 


K  =  2  _ _ 

L  I  ST  ="”o 

L  INDA  __=  1  _  _ _ 

^  D  =  Win  -  W'lN) 

00  55_  I^tN  _  _  _ 

W(I)  =  100.-(100.*(W( I )-W(N) )/D) 

c  0  N  V  E  rT  tip  e”1  n  Tn  c  hFs  'To  minutes^ 

■■  '"ti  I )  =rTTl  )  /sp£'ed 
5  5  ,_a3NTiNyE _ _ _ _ 

ID T  =_  P  E  N  S  £  PARA  T  1 0  N  I _N  ._M_I  NU  T  E  S^__. . 

0  T  _  =  ,  S  P  E_EJD  *  0  F  F  S  E_T _ _ _ _ 


..STABI,  NA_JQP  .L00P._. .  _ _ _ 

_D0  ICO  NW  =^,99  _  _ 

5"8  II  =  LiNDA-l 

_  _  WW(NW)  =  FLOAT  )  _ _ _ _ _ 

SCAN  WEIGHT  DATAFGR  ONE  CLOSE  TO  BUT  JUST  GREATER  THAN  ONE  PERCENT  WEIGHT _ 

LOSS.  II  =rN0E'X  OF  THAT  POTnT 

D(D^60“  r=rrN  DA  ,lN  . 

11=11+1  _  _  _ _ _ 

"  IF  (W  (  I  ).Gt.WW"(Nw')  )  G(D  TO  70 

60  CONLINUE  _  ___  _ 

70  LINDA  =  II-(LL/2) 

lT'n'DA” ="' ‘i  ND E X  0 F~F  I R S f  '  0 ATA  TO  BE  USTD  BY  PLSQ 
bo’  6c"j=i7lT  . 

_JI  =  LINDA+J-l  _ 

X  (  J  )  =  T  (  J  I  ) 

_  Y  (  J  )  =  W  (  J  I  )_ _ 

Z"(  J)  =  TEMPI  JI  ) 

lC_CJD]iTJ  NJJE _ _ 

_ CURVE  FIT  CF  TIME  AND  WEIGHT  DATA _ _ 

_ CALL  PLSQ  (X,Y,NN,K,C,LIST,EMAX,ERMS,EMEQ) _ LM 

KK  =  1 

STT^T"  TbOP”'’TO^ CHECK  FOR  BAD  INPUT  DATA 
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C 

DO  61  J=1,LL 
JI  =  LINDA+J-1 


C 

C 

WE  =  WEIGHT  CALCULATED  FROM  POLYNOMIAL 

c 

WE  =  C(  1)*T(  jn**2+C(2)*T(  JI  )+C(3) 

c 

_ ^ 

COMPARE  CALCULATED  AND  ORIGINAL  DATA 

c 

IF  {ABS(WE-W(JI)).GT.l.)  GO  TO  82 

82 

U  1  L  b  1 

WRITE  (6,4000)  JI,W(jn,WE 

c 

c 

REPLACE  BAD  DATA  BY  CALCULATED  VALUES 

c 

W(JI)  =  WE 

81 

KK  =  2  — - 

CONTINUE 

c 

GO  TC  (83,58  ),KK  "  '  . .  . 

c 

c 

CHECK  FOR  IMAGINARY  ROOTS  IN  SOLUTION  OF  QUADRATIC 

83 

SCREW  =  C(2)*C(2)-4.0*C(1)*(C(3»-WW(NW))  . . . .  -  - 

IF  (SCREW.lt. 0.0)  GO  TC  90 

c 

c 

r 

USE  REAL  RCOT  TO  DETERMINE  TIME  CORRESPONDING  TO  EACH  PERCENT  WEIGHT  LOSS 

TNW(NW)  =  (SQRT(C(2)*C(2)-4.0»C(1)*{C(3 >-WW(NW ) ) )-C ( 2 ) ) /( 2. 0*C ( 1 ) ) 


C  OWDT  =  RATE  OF  WEIGHT  LOSS 


DWOT(NW)  =  2.0*C(1  )nNW(NW)  +  C(2) 

126 

GO  TC  100 

. 

90  TNW(NW)  =  T(II-l) 

DWDT(NW)  =  0.0 

TA  =  T(II)*SPEED 


C 

C 

WRITE  OUT  IDENTIFICATION  AND  LOCATION 

OF  BAD  DATA 

c 

WRITE  (6,3160) 

NW,II,TA,W(II) 

136 

100 

CONTI NUE 

K  =  5 

c 

c 

CURVE  FIT  CF  TIME  AND  TEMPERATURE  DATA 

c 

CALL  PLSQ(T,TEMP,N,K,C,LIST,EMAX,ERMS, 

EMEQ) 

142 

1 - 

WRITE  (6,3100) 
WRITE  (6,5200) 

EMAX 

ERMS 

143 

144 

( 

WRITE  (6,5300) 
WRITE  (6,5400) 

(C(I),I=1,6) 

145 

146 

■ 

WRITE  (6,3110) 
DO  120  NW=1,99 

153 

c 

c 

CT  =  CORRECTED 

TIME,  USED  TO  FIND  TIME 

AT  EACH  TEMPERATURE 
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CT  =  TNW(NV<)-DT 

C 

TSTOR  =  C(l) 

C 

C 

LOOP  TO  EVALUATE  TEMPERATURE  POLYNOMIAL  FOR  EACH  VALUE  OF  CT 

200 

DO  200  1=2,6 

TSTOR  =  TSTOR*CT+C( I ) 

TPOLY(NW)  =  TSTOR 

TSTOR  =  5. *0(1) 

DO  250  1=2,5 

J  =  6-1 

C 

250 

TSTOR  =  TSTOR*CT+FLOAT{ J)*C( n 

C 

C 

TDER  =  TEMPERATURE  DERIVATIVE 

RTEMP  =  RECIPROCAL  OF  ABSOLUTE  TEMPERATURE 

C 

TDER(NW)  =  TSTOR 

RTEMP(NW)  =  1.0/(TPOLY(NW)+273.16) 

WRITE  (6,3120)  NW, OWDT (NW ), TPOLY ( NW ) ,TDER ( NW ), RTEMP (NW ) 

172 

120 

CONTINUE 

STDER  =  0.0 

C 

C 

CALCULATE  AVERAGE  TEMPERATURE  DERIVATIVE  (AVE) 

C 

DO  125  1=1,99 

125 

STDER  =  STDER  +  TDER(I) 

CONTINUE 

AVE  s  STDER/99.0 

WRITE  (6,3125»  AVE 

C 

C 

SET  UP  DUMMY  POINTS  FOR  GRAPH  PLOTTING  SUBROUTINE  ( GP ) 

C 

WW(IOO)  =  0.0 

187 

DWDT(IOO)  =  0.0 
'  TNW(IOO)  =  TNW(99) 

TDERllOO)  =  0.0 

WH(lOl)  =  100.0 

DWDT(lOl)  =  0.0 

TNW(lOl)  =  TNW(99» 

TDER(lOl)  =  TDER(99) 

WRITE  (6,3000) 

188 

WRITE  (6,3130)  ID 

L  =  3 

189 

LS  =  5 

LW  =  101 

LN  =  50 

M  =  101 

DATA  A/IH./ 

JN  =  1 

C 

C 

PLOT  GRAPH  OF  RATE  OF  WEIGHT  LOSS  AGAINST  PERCENT  WEIGHT  LOSS 

C 

CALL  GP  (WW,OWDT,L,LS,M, JN,LW,LN,A^PLOT) 

196 

WRITE  (6,3000) 

WRITE  (6,3140)  ID 

197 

C 
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c 

PLOT  GRAPH  OF  PERCENT  WEIGHT  LOSS  AGAINST  TIME 

c 

CALL  GP  (TNW,WW,L,LS,M,JN,LW,LN, A.PLOT ) 

198 

199 

WRITE  (6,3000) 

WRITE  (6,3150)  ID 

200 

c 

c 

PLOT  GRAPH  OF  TEMPERATURE  DERIVATIVE  AGAINST  TIME 

c 

CALL  GP  (TNW,TDER,L,LS,M, JN,LW,LN,A,PLOT) 

201 

202 

DWDT(IOO)  =  0.0 

TPOLY(IOO)  =  0.0 

C 

c 

PUNCH  OUTPUT  CARDS  CONTAINING  PERCENT  WT.  LOSS(NW)  THEN  THREE 

PAIRS 

OF 

c 

c 

temperature  and  rate  of  weight  LOSS  DATA 

DO  150  NW=l,100t3 

PUNCH  5000,ID,NW,DWDT(NW) ,TPaLY(NW),DWOT(NW+l) ,TP0LY(NW+1)  , 

150 

.DWDT(NW+2) ,TP0LY(NW+2) 

CONTINUE 

206 

C 

C 

LOOK  FOR  FURTHER  SETS  OF  DATA  IF  NEXT  CARD  CONTAINS  A  ONE  IN 

COLUMN 

10 

C 

READ  (5,1020)  MORE 

215 

IF  (KORE.EC.l)  GO  TO  1 

STOP 

1000 

i 

F0RMAT(8X, A4,2X,A6,A2,2X, A6, A6, A6,A3,F7.4,2X,F5.4,2X,F6.2,3X, 
,I3,3X,I3) 

1010 

1020 

F0RMAT(13X,I1,2X,3(F6.2,F5.1,1X,F4.0),F6.2,F5. 1,1X,F4.0) 
FORMAT  (9X,I1) 

1030 

,3000 

FORMAT  (12) 

FORMAT  (IHl) 

3010 

3050 

FORMAT  (10X,21HT0TAL  NO  OF  POINTS  =  ,14) 

F0RMAT(5X, A4,8X,A6,A2, lOX, A6, A6, A6,  A3//) 

3060 

3065 

FORMAT  (10X,33HCHART  SPEED  (INCHES  PER  MINUTE)  =,F6.4) 

FORMAT  ( 10X,38HCHART  SPEED  (FURLONGS  PER  FORTNIGHT)  =,F7.4) 

3070 

3080 

FORMAT  (10X,17H0FFSET  (INCHES)  =,F7.4) 

FORMAT  (10X,17HXRANGE  (INCHES)  =,F7.2) 

3090 

3100 

FORMAT  (10X,18HINITIAL  TEMP  (C)  =,I5) 

FORMAT  (10X,16HFINAL  TEMP  (C)  =,I5) 

3106 

3110 

FORMAT  (18X, I5,4X,F9.3,4X,F7.2,4X,F8.2) 

FORMAT  (//3X,11HWEIGHT  LOSS , 6X , 8HDWDT ( NW ) , 1 3X , 5HTP0L Y , 6X, 

3120 

,4HTDER,11X,5HRTEMP) 

FORMAT  (6X,I3,10X,E12.5,7X,F9.3,2E15.5) 

3125 

3130 

FORMAT  (//10X,27H  AVERAGE  TEMP  DERIVATIVE  =  ,E15.5) 

FORMAT  (10X,19H0WDT  VS  WEIGHT  L0SS,20X,A4) 

3140 

3150 

FORMAT  (10X,19HWEIGHT  LOSS  VS  TIME,20X,A4) 

FORMAT  (10X,12HTDER  VS  TIME,20X,A4) 

3160 

F0RMAT(2X, 17HSCREW  LESS  THAN  0 , 1 OX , 3HNW= , I 3 , lOX , 3HI I = , I 4, lOX, 
2HT=,F6.2,10X,2HW=,F5.1) 

3170 

4000 

FORMAT  (10X,25HN0  OF  PTS  IN  CURVE  FIT  =  ,12) 

FORMAT  (10X,9HAT  PT  NO  ,I4,10H  WEIGHT  =  ,F5.1,13H  REPLACED  BY 

f 

• 

5000 

F5.1)  ^ 

FORMAT  (2X,A4,I4,E13.5,F6.1,E13.5,F6.1,E13.5,F6.1) 

5100 

5200 

FORMAT  (10X,17HMAX  TEMP  ERROR  =  ,F10.6) 

FORMAT  (10X,30HTEMP  ROOT  MEAN  SQUARE  ERROR  =  ,F10.6) 

5300 

FORMAT  (10X,15HTEMP  POLY  COEFF) 

TGA 

PLOT  -  EFN  SOURCE  STATEMENT  -  IFNIS)  - 

5400  F0RMAT(13X,FI0.6) 
END 
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'’LOT  -  EPN  SOURCE  STATEMFiviT  -  IFN(S)  - 

C 

in  '16  1  iNt-  IGA  Ra^AMsTERS  BY  F  R  1  E  0  MANTSMFT  HTO' 

programme  accepts  data  capos  having  three  sets  of  OATA  P'^R  card. 

c 

c 

LAST  CARD  OF  rACH  DECK  MUST  HAVE  A  OMt  IN  COLUMNl.  LAST  CARD  OF 

LAST  DECK  FOR  ONE  pOLVMFR  SYSTEM  MUST  HAVE  A  TWO  IN  COLUMN  1  INSTEAD 

u 

c 

TO  RUN  A  ScCOND  SFT  OF  DECKS.  PUNCH  A  CARD  WITH  A  THR^E  IN  COLUMN  1 

AND  PLACE  between  SETS 

0 

c 

AT  THE  TND  (TF^'ALL  DECKS  PLATE  A  BLANK  CARD  THEN  ATI  FFDF - - - - 

c 

c 

SYMBOLS  DWO"^  -  RATE  OF  WEIGHT  LOSS,  RTEMd  =  RECIPROCAL  OF  ABSOLUTE 
TEMPFRATUPF,  RATE  =  LOG  RATE  OF  WEIGHT  LOSS,  SLOPE  =  SLOPE  OF  ARRHENIUS 

c 

c 

c 

c 

plot,  prfx  -pRE-EXPONsNI  IAL  factor,  plot  -  DIM'-NSII'N  OF  GP  SUBRO'UTINF  . 

~  ACTIVATION  =NEPGY,  X  AND  Y  REPRESENT  DATA  TREATED  BY  GP 

!  f^ULY  INPUT  r  ATUR  „S  f  10  -  I  I  ^  T  C  AT  I  ON ,  A  =  NO,  OF  SYM5DU>  P^T  fTP — 

AA  =  PCPCENT  HEIGHT  LHSSt  ^'FW  =  FUMCTIOM  FPOM  FPIEO^ANS  EQUATION 

c 

c 

FW  AVFRAGE  AFW,  bB  -  LOG  (PERCFNT  PFSIDUF),  WF  =  AVFRAGE  AFW  ‘ 

DIMENSION  OWDT ( 1 OC , 1C ) , R  TE mp ( 1 QO , 1 0 ) , R A TF ( 1 00 , 1 0 )  .SLOPF ( ICO ) , 

.PREXnOO)  .PLOT  (50, 100)  ,  ACTE  (  1 00 )  ,  X  {  10  )  ,  Y  (  1 0 )  , 

•  f  P0LV(  100, 10) ,  Tn(  10 ) , A ( 1 ) , AA( 100)  , AFW( 1 0) ,FW( 100 ) , BR( =5 ) ,WF ( 05 ) , 
.SPS(100),SDS( IOC) ,SDI( 100) 

1 

read  (5,1C0C)  IG, COMl, COM2, CQM3, COMA, C0M5, C0M6  ,C0M7, CMS 

WRITE  (6,3000) 

1 

3 

? 

10 

write  (6,1100)  10, COM!, COM2, COMO, COMA, CDM5,CnM6, COM?, CCM8  ■■ 

J  =  0 

4 

c_ 

J  =  J  +  1  . . - 

c 

c  _ 

START  LOOP  TO  REAC  IN  OATA  - - - - — 

c 

c 

c__ 

DO  20  NW  =  1,  ?■',?  - - - 

LBj  =  1  IN  COLUMN  1  OF  LAST  CARD  OF  A  DECK,  LAST  CARD  OF  LAST  Dp'CK  FDR 

ONE  POLYMER  SYSTEM  NEEDS  LBJ  =  2. 

c 

RFAD  (f,  1200)  LBJ,  in(  J  ),  IW,OWOT(NW,J),TPnLY(NW,J)  ,  DW-ST  (  NW+ 1 ,  J  )  , 

c 

c 

c 

•  TPOLYl  NW+!  ,  J  )  ,  D'WDTI  NW+2,  J  )  ,  TPOLY(  NW  +  2,  J)  ' - 

cHbtK  1  HA  i  iN^UI  L A<-iL)S  A^E  IN  CONSECUTIVE  npoF^*'  *'  *"  - - 

10 

3 

1  H  (1  W—  M  W  )  •?  f  4  f  1  — - — 

WRITE  (6,19^0)  NW ,  I 0 ( J  )  , I W 

21 

4 

STOP  - - - - 

AA(NW)  =  FLOAt(nW) 

AA(NW+?)  =  FL0AT(NW+2) 

?0 

IF  (LBJ.EQ.l)  GO  Tn  ir  - — - 

I^(L3J.F0.2)  GO  TO  25 

c 

c 

c_ 

2  5 

X  J  —  J 

write  list  of  run  IDS 

r 

WRITF  (6,  !  800  )  (  ID(  I  )  ,  1  =  1,  j1  - - - - - 

c 

c 

che{:k  for  at  L'-as"^  three  data  dfcks 

T  4 

t  R  (  J-  3  )  3  C  ,  3  5 , 3  5  - 
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30 

WRITE  (6,2000) 

44 

35 

GO  TO  300 

WRITF  (6.1509) 

46 

TSU.M  =  0 

N  =  C 

C 

SPREX  =  0.0 

c 

c 

START  LOOP  TO  CALCULATE  LEAST  SQUARES  LIME  OF 

LOG (PATE)  VS. 

PTtMP 

DO  45  MW  =  4,08 

SUMXX  =  0  .  . 

SUMYY  =  0 

SUMX  =  C 

SUMY  =  0 

SUMXY  =  0  _ _ - 

C 

DO  40  K  =  1,J 

C 

c 

CHECK  FOR  7EP0  RATES 

lT(bW0T(N’w,K)  .LT.l.OE-lO)  GO  TO  65 

RATF(MW,K)  =  ALOG10(DWDT(MW,K) ) 

53 

RTEMP(MW,K)  =  1.0/(TPnLY(NW,K)+273.16) 
C _ 


SUMXX  =  05STI4L  SUM  OF  X  SOUAPED  ETC- 

“"“SUm'xx'  =''SUMXX  rTf=MprNw7KT*^2  . . 

_ SUMYY  =  SUMY Y  (  P  A  TF  (  NW  ,  K  )  )  **  2 _ 

SUMX  =  SUMX  +  RTEMP(NW,K) 

_  SUMY  =  SUMY  *  RATF(NW,K) _ 

49  SUMXY  =  SUMXY  +  RTEMP(NW,)<)’!'RATF(  NW,K) 

_f;r)  TO  5  5  _ _ _ _ _ _ _ 

C . .  . 

C  _ SET  UP  oijyyy  PHINTS  FOP  9P  IF  A  UWDT  VALU^  IS  ?FPP 

C 

3 5  ACTE(N!W)  =  0. _ 

PRFX(NW)  =  b. 

CATE  (NW,K)  =  0. _ _ _ _ _ _ _ . 

. . '"rTFMP  (NW,K)  =  0.0015 


GO  TO  45  _ _ _ _ 

55  SLOPr(NW)  =  (XJ*SUMXY-SUMX=»SUMY)/(  XJ7.SUMXX-SUMX^'*2) 

S  PS  (  \iw  )  =(  (  SUMYY-  (  SUMV«SUMY/X.I  )-•  (  (  XJ* SUMXY-  S')'-* X* SU^Y )  / 

.(  XJ*XJ*SUMXX-XJ*SUMX*SUMX)  )  ) /( XJ-2.C)  ) 

ALPHA  =  (SPS(NW)/(S  IJMXX- (  SUMX*SUMX/XJ  )  )  )  *4. 576 _ 

I F ( A  L  P H A )  5  8 , 5  8 , 57 

57  SDS(NV))  =  SORT  (ALPHA) _ 

GO  TO  50 

33  SDS(NW)  =  9.0 _ 

59  beta  =  (  SPS  (  MW  )  *SUMXX/ (  XJ*  SUMX  X- SUMX^SU'^X  )  ) 

_IJiiRFTA  )  62t  62t  61 _ _ _ 

61  S  0  I ( MW j  =  SOPT(RETA) 

GO  TO  63 _ 

62  SOI  (MW)  =  0.0 

63  ACTF(NW)  =  - SL QP E ( MW ) * 4 . 57  6 _ 

PPFX (MW  )  =  ( SUMXX* SUMY- SUMX*SUMXY ) / ( X J*S UMX X- SUM X ) 

LP  (JMW.LT^10J_Gn_TO  « _ _ 

IF(NW.GT.BC)  GO  TP  45 
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TSUM  =  TSUW-SLOPE(NW) 
SPREX  =  SPftEX  +  Pfiex(Ni 
N  =  N+1 
“CONTINUE 


SOURCE  STATEMENT  -  IFNIS) 


AVPREX  =  SoREX  /  FLOAT(N) - 

AVEA  =  TSUM/FL0AT(N) 

AVACTE  =  AVEA*1. 987*2. 303  “ 

START  LOOP  TO  CALCULATE  AFW 

“DO  70  NW  =  4, P'S  — 

Z  =  0 

00  <50  K  =  1,J  ' - 

AFW(K)  =  RATE{NW,K)  +  A VEA^RTP ( NH t K J 

Z  =  Z  +  AFW(K)  - 

FW(NW)  =  Z/XJ 

WN  =  FLOATINW)  “ - 

C,G  =  ALOG10(  lOO.-WN  ) 

SD  =  0  - 

DO  93  K  =  1,J 

SO  =  SO  (FW(NW)~AFW(K) 


SQRTISD/YK) 


c 

_ ^ 

WRITE  OU^  RESULTS  PERCENT  WT.  LOSSi  ACTIVATION  RNFPGY,  PPE-FXPONENT lAL 
FACTOR,  AVERAGE  FW,  AND  STANDARD  DEVIATIONS,  ALSO  LOG  WEIGHT  PFMA I N TNG{ GG » 

c 

70  WRITE  (6,1400)  *^W,  ACTE  (  NW )  ,  SDS  (NW  )  ,  PR  E  X  (NW )  ,  SD  I  (  NW  )  ,  FW  (  NW )  .  SDAFW  , 

13  3 

.GG  ■"  ■"  -  - — - 

WRITE  (6,1425)  AVACTE 

134 
14  1 

WKlIt{6, 1435)  AVPREX  '  - — 

WRITE  (6,1440) 

14  2 

c 

_ c 

SET  UP  INFORMATION  FOR  Gp  SUBROUTINE,  SEE  OTHER  PROGRAMS 

c 

L  =  3 

LS  -  5  - - - 

LW  =  100 

L  N  —  50  - — - - 

M  =  J 

0  A  T  A  A  /  i  H .  /  .  ‘  -  - - - 

JN  =  1 

c 

c 

START  LOOP  FOR  PLOTTING  GRAPHS  AT  10  PFRC^NT  WEIGHT  LOES  INTprvaIs 

c 

00  200  NW  =  10,99,10 

_ 

DO  100  K  -  1,J  - 

X(i<)  =  RTEMP(NW,K) 

iuu  YIM  -  KAIb(NWfK)  - - - - 

WRITE  (6,3000) 

16C 

c 

WRITE  (6,1700)  NW  . . - 

c 

PLOT  GRAPH  OF  LOG  (RATE  OF  WEIGHT  LOSS)  AGAINST  PECIPPGCAL  '  “ 
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C 

OF  TFMPFFATUPF 

C 

300 

CALL  GP  (X,Y,L,LS,M,JN,LW,LN,A,PLOT) 

16  1  ' 
163 

M  =  100 

WRITE  (6,30C0) 

167 

C 

WRITE  (6,3100) 

c 

c 

PLOT  GRAPH  OF  ACTIVATION  ENERGY  AGAINST  PERCENT  WT.  LOSS 

168 

CALL  GP  (AA,  ACTE,L,LS,M, JN,LW,LN,A,PL0T) 

WRITE  (6,3000) 

15S 

17C 

c 

WRITE  (6,3200) 

c 

c 

PLOT  GRAPH  OF  po E~ EXPONENT  I AL  FACTOR  AGAINST  PERCENT  WEIGHT. LOSS 

171 

CALL  GP  (  AA,prEX,L,LS,M, JN,LW,LN,A,PLOT) 

WRITP  (6,3000) 

172 

173 

WRITE  (6,3300) 

DO  75  1=1,87 

174 

7F 

88(1)  =  ALaG10(10C.-AA( 1+3)) 

WF{  I  )  =  FW( 1  +  3 ) 

179 

LW  =  P5 

M  =  87 

C 

c 

PLOT  GRAPH  OF  LOG(AFW)  AGAINST  LQG(pEFCENT  RESIDUE  WEIGHT) 

c 

CALL  GP(BB,WF,L ,LS, M, JN,LW,LN, A,PLOT) 

188 

c 

WRITE  (6,3000) 

c 

c 

LOOK  F<PR  FURTHER-  SETS  OF  DATA 

189 

300 

READ  (5,1300)  MORE 

IF{mgrf.=q.3)  go  TO  1 

19C 

PO 

1000 

STOP 

format  ( 2X , A3, 2X, eA6) 

1100 

1200 

FORMAT  (  lOX, A3, 2X, 8A6) 

FORMAT  (I1,1X,A4,I4,E13.5,F6.1,P13.5,F6.1,F13.5,F6.T) 

1300 

1400 

FORMAT  (ID 

FORMAT  ( lOX, n ,4X,-3PF7. 3, 5X,F6. 3, 5X,CPF5.3,2 ( 5X,F6.3)  ,2 ( 5X,F6.4) ) 

1425 
143  5 

FORMAT  (//10X,39H  average  ACTIVATION  ENERGY  =  ,-3PF6.3) 

FORMAT  (1CX,17H  AVERAGE  LOG  PREX,10X,2H=  ,F6.3) 

1440 

1500 

FORMAT! lOX, 34H80TH  FOR  10-80  PERCENT  WEIGHT  LOSS) 

FORMAT  (/8X,7HWT  LO S S , 2 X , 8  HE A ( KC A L ) , 3X , 8H ST . OE VN . , 3 X , 8H LOG  PREX, 

1700 

.3  X, BHST.OEVN. , 2X, 10HAV.LOG  AFW , 2X , BHST . DE VM . ,2X,11HL0G  RES.WT.  ) 

FORMAT  (10X,18HLnG  RAT'^  VS  1 /TEMP  /  lOX  ,  14HWE  I GHT  LOSS  =  ,14) 

1800 

1900 

FORMAT  (/10X,11HRUN  ID  NOS  ,9(A4,2H,  )) 

FORMAT  ( ICX, 13HEPR0R  FOR  W  =,I4,7HRUN  NO  ,A3,eH  READ  ,13, 

2000 

.9H  INSTEAD.) 

FORMAT  (1CX,25HLFSS  THAN  3  HEATING  RATES/lHl ) 

300  0 
310  0 

FORMAT  (IHl) 

format  ( 1GX,32HACTIVAT ion  ENERGY  VS  WEIGHT  LOSS) 

37CO  FORMiT  (  10'X,?2HPR'^~FXP  VS  WEIGHT  LOSS) 

3300  FnP«^4T(  10X,46HAVFP  LOG  AF(W)  VS  LOG  PERCENT  WEIGHT  REMAINING) 


FNID 
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ICTF  19/03/68  TEFLON  (NEW  BAL.) 

CHART  SPEEC  (INCHES  PER  MINUTE)  =1.0000 
CHART  SPEEC  (FURLONGS  PER  FORTNIGHT)  =  2.5455 
OFFSET  ( INCHES)  =  0.0625 
XRANGE  (INCHES)  =  36.60 

INITIAL  TEMP  (C)  =  418 

final  TEMP  (C)  =  625 

NO  OF  PTS  IN  CURVE  FIT  =  10 
TOTAL  NO  OF  POINTS  =  200 

MAX  TEMP  ERROR  =  2.616966 

TEMP  ROOT  MEAN  SQUARE  ERROR  =  0.635928 

TEMP  POLY  CDEFF 
0.000003 
-0.000230 
0.007066 
-0.087217 
5.973110 
418.330875 


WEIGHT  LOSS 
1 
2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 


CWOT(NW) 
0.43704E-00 
0.76467E  00 
0.10531E  01 
0.13932E  01 
0.17160E  01 
C.19640E  01 
0.22448E  01 
0.25599E  01 
0.29225E  01 
C.32235E  01 
0.34049E  01 
0.35286E  01 
C.37758E  01 
0.40299E  01 
0.42162E  01 
0.44440E  01 
0.46574E  01 
0.50011E  01 
0.53156E  01 
0.55268E  01 
0.56799E  01 
0.57782E  01 
0.58077E  01 
0.59939E  01 
0.61520E  01 
0.63133E  01 
0.65319E  01 
0.68436E  01 
0.70771E  01 
0.73456E  01 
0.74898E  01 
0.75789E  01 
0.76239E  01 
0.76522E  01 
C.77316E  01 
0.78633E  01 
0.80040E  01 
C.82437E  01 
0.83949E  01 


TPOLY 

507.899 

518.102 

524.516 

529.237 

532.890 
535.986 

538.693 
541.067 
543.174 
544.993 
546.670 
548.339 

549.891 
551.333 

552.693 
553.998 

555.237 
556.430 

557.516 
558.527 
559.508 
560.473 
561.465 
562.406 
563.326 
564.223 
565.113 
565.958 
566.750 
567.515 
568.251 
568.978 
569.701 
570.446 
571.172 
571.894 
572.598 
573.284 
573.950 


TDER 

0.57C00E  01 
0.57125E  01 
0.57103E  01 
0.57C40E  01 
0.56964E  01 
0.56e84E  01 
0.56803E  01 
0.56724E  01 
0.56649E  01 
0.56581E  01 
0.56516E  01 
0.56449E  01 
0.56385E  01 
0.56326E  01 
0i56269E  01 
0.56214E  01 
0.56162E  01 
0.56113E  01 
0.56C68E  01 
0.56C27E  01 
0.55988E  01 
0.55950E  01 
0.55912E  01 
0.55877E  01 
0.55843E  01 
0.55812E  01 
0.55782E  01 
0.55754E  01 
0.55729E  01 
0.55707E  01 
0.55686E  01 
0.55666E  01 
0.55648E  01 
0.55631E  01 
0.55615E  01 
0.55601E  01 
0.55588E  01 
0.55577E  01 
0.55568E  01 


RTEMP 

0.12803E-02 

0.12638E-02 

0.12536E-02 

0.12463E-02 

0.12406E-02 

0.12359E-02 

0.12318E-02 

0.12282E-02 

0.12250E-02 

0.12223E-02 

0.12198E-02 

0.12173E-02 

0.12150E-02 

0.12129E-02 

0.12109E-02 

0.12090E-02 

0.12072E-02 

0.12054E-02 

0.12038E-02 

0.12024E-02 

0.12010E-02 

0.11996E-02 

0.11981E-02 

0.11968E-02 

0.11955E-02 

0.11942E-02 

0.11929E-02 

0.11917E-02 

0.11906E-02 

0.11895E-02 

0.11885E-02 

0.11875E-02 

0.11864E-02 

0.11854E-02 

0.11844E-02 

0.11834E-02 

0.11824E-02 

0.11814E-02 

0.11805E-02 
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40 

0.85006E 

01 

41 

0.85905E 

01 

42 

C.87186E 

01 

43 

C.87468E 

01 

44 

0.87688E 

01 

45 

0.88183E 

01 

46 

0.88743E 

01 

47 

0.91338E 

01 

48 

0.92862E 

01 

49 

C.94210E 

01 

50 

C.95533E 

01 

51 

0.97158E 

01 

52 

C.97623E 

01 

53 

0.97623E 

01 

54 

0.97445E 

01 

55 

0.97178E 

01 

56 

C.99691E 

01 

57 

C. 10038E 

02 

58 

0. 10131E 

02 

59 

C. 10208E 

02 

60 

0.  10281E 

02 

61 

C. 10400E 

02 

62 

C.10414E 

02 

63 

C. 10376E 

02 

64 

C.  10310E 

02 

65 

0. 10204E 

02 

66 

C.10352E 

02 

67 

0. 10344E 

02 

68 

0. 10358E 

02 

69 

0. 10391E 

02 

70 

0.10387E 

02 

71 

C. 10423E 

02 

72 

0. 10433E 

02 

73 

0. 10321E 

02 

74 

0.  10237E 

02 

75 

C.  10099E 

02 

76 

0. lOOllE 

02 

77 

0. 10151E 

02 

78 

0. 10124E 

02 

79 

0. 10033E 

02 

80 

C.99414E 

01 

81 

0.97937E 

01 

82 

0.96724E 

01 

83 

C.94010E 

01 

84 

C.90958E 

01 

85 

0.86745E 

01 

86 

0.84766E 

01 

87 

0.82320E 

01 

88 

0.79915E 

01 

89 

0.77573E 

01 

90 

0.75096E 

01 

91 

0.72152E 

01 

92 

C.68785E 

01 

93 

0.64873E 

01 

94 

C.59624E 

01 

95 

0.54449E 

01 

96 

C.48352E 

01 

97 

0.41379E 

01 

98 

0.33280E 

01 

99 

0.23131E 

01 

AVERAGE  TEMP  DERIVATIVE  = 


574.597 

0.55560E 

01 

0.11796E-02 

575.245 

0.55554E 

01 

0.11787E-02 

575.883 

0.55549E 

01 

0.11778E-02 

576.515 

0.55545E 

01 

0.11769E-02 

577.151 

0.55543E 

01 

0.11760E-02 

577.783 

0.55543E 

01 

0.11752E-02 

578.419 

0.55543E 

01 

0.11743E-02 

579.047 

0.55546E 

01 

0.11734E-02 

579.648 

0.55550E 

01 

0.11726E-02 

580.236 

0.55555E 

01 

0.11718E-02 

580.817 

0.55562E 

01 

0.11710E-02 

581.385 

0.55570E 

01 

0.11702E-02 

581.946 

0.55579E 

01 

0.11694E-02 

582.510 

0.55590E 

01 

0.11687E-02 

583.086 

0.55602E 

01 

0.11679E-02 

583.674 

0.55617E 

01 

0.11671E-02 

584.248 

0.55632E 

01 

0.11663E-02 

584.802 

0.556495 

01 

0.11656E-02 

585.351 

0.55668E 

01 

0.11648E-02 

585.896 

0.55687E 

01 

0.11641E-02 

506.432 

0.55708E 

01 

0.11633E-02 

586.964 

0.55731E 

01 

0.11626E-02 

587.493 

0.55755E 

01 

0.11619E-02 

588.031 

0.55780E 

01 

0.11612E-02 

588.572 

0. 558085 

01 

0.11605E-02 

509.117 

O.55038E 

01 

0.11597E-02 

589.667 

0.558705 

01 

0.11590E-02 

590.210 

0.55903E 

01 

0.11583E-02 

590.751 

0.55938E 

01 

0.11575E-02 

591.290 

0.55975E 

01 

0.11568E-02 

591.825 

0.56C13E 

01 

0.11561E-02 

592.364 

0.56C54E 

01 

0.11554E-02 

592.894 

0.56C96E 

01 

0.11547E-02 

593.423 

0.56140E 

01 

0.11540E-02 

593.962 

0.56186E 

01 

0.11532E-02 

594.509 

0.56235E 

01 

0.11525E-02 

595.089 

0.56290E 

01 

0.11517E-02 

595.659 

0.56346E 

01 

0.11510E-02 

596.219 

0.56403E 

01 

0.11502E-02 

596.774 

0.56463E 

01 

0.11495E-02 

597.328 

0.56524E 

01 

0.11488E-02 

597.889 

0.56589E 

01 

0.11480E-02 

598.456 

0.56657E 

01 

0. 11473E-02 

599.041 

0.56729E 

01 

0.11465E-02 

599.649 

0.56808E 

01 

0.11457E-02 

600.311 

0. 56897E 

01 

0.11449E-02 

600.964 

0.56988E 

01 

0.11440E-02 

601.653 

0.57C88E 

01 

0.11431E-02 

602.368 

0.57196E 

01 

0.11422E-02 

603.100 

0.57311E 

01 

0.11412E-02 

603.846 

0.57433E 

01 

0.11402E-02 

604.618 

0.57565E 

01 

0.11392E-02 

605.430 

0.57709E 

01 

0.11382E-02 

606.282 

0.57866E 

01 

0.11371E-02 

607.219 

0.58C47E 

01 

0.11359E-02 

608.232 

0.58252E 

01 

0.11346E-02 

609.365 

0.58493E 

01 

0.11331E-02 

610.678 

0.58787E 

01 

0.11314E-02 

612.262 

0.59164E 

01 

0.11294E-02 

614.390 

0.59709E 

01 

0.11267E-02 

0.56306E  01 
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TEFLON  DEGRADATION 

H.R.  =  HEATING  RATE  IN  DEGREES  C  PER  HOUR 
RATE  CF  WT.  LOSS  IN  PERCENT  PER  MINUTE 


I  .D./H.R. 

4TF/ 

85 

9TF/145 

llTF/44 

12TF/ 

55 

PERCENT 

RATE  OF 

TEMP. 

RATE  OF 

TEMP. 

RATE  OF 

TEMP. 

RATE  OF 

TEMP. 

WT.LOSS 

WT.LOSS 

DEG  C 

WT.LOSS 

DEG  C 

WT.LOSS 

DEG  C 

WT.LOSS 

DEG  C 

1 

0.074 

486.5 

0.207 

494.9 

C.069 

480.8 

0.078 

479.9 

2 

0.210 

497.6 

0.390 

505.3 

0.138 

489.4 

0.154 

489.4 

3 

0.312 

503.4 

0.553 

511.2 

0.194 

494.5 

0.212 

494.6 

4 

0.418 

507.6 

0.708 

515.6 

0.243 

498.3 

0.265 

498.5 

5 

0.526 

510.8 

0.857 

519.0 

C.307 

501.3 

0.323 

501.7 

6 

0.595 

513.4 

0.998 

521.9 

0.349 

503.6 

0.383 

504.3 

7 

0.692 

515.8 

1.134 

524.4 

C.403 

505.8 

0.433 

506.6 

8 

0.776 

517.8 

1.286 

526.6 

0.451 

507.6 

0.504 

508.7 

9 

0.825 

519.6 

1.412 

528.5 

C.509 

509.2 

0.594 

510.4 

10 

0.910 

521.3 

1.590 

530.2 

0.563 

510.7 

0.634 

511.9 

11 

0.987 

522.8 

1.722 

531.8 

0.598 

512.0 

0.637 

513.4 

12 

1.091 

524.3 

1.811 

533.2 

C.640 

513.3 

0.688 

514.8 

13 

1.185 

525.5 

1.932 

534.6 

C.688 

514.4 

0.732 

516.1 

14 

1.239 

526.7 

2.054 

535.8 

0.720 

515.5 

0.789 

517.3 

15 

1.285 

527.8 

2.193 

537.0 

0.773 

516.5 

0.838 

518.4 

16 

1.363 

529.0 

2.327 

538.1 

0.805 

517.5 

0.  882 

519.5 

17 

1  .451 

530.0 

2.398 

539.2 

0.840 

518.4 

0.921 

520.6 

13 

1.532 

530.9 

2.513 

540.3 

0.868 

519.3 

0.981 

521.6 

19 

1.622 

531.9 

2.678 

541.2 

0.948 

520.1 

1.028 

522.5 

20 

1.704 

532.7 

2.813 

542.1 

0.989 

520.9 

1.066 

523.4 

21 

1.738 

533.6 

2.873 

543.0 

1.009 

521.6 

1.107 

524.2 

22 

1.793 

534.4 

2.878 

543.9 

1.023 

522.3 

1.146 

525.1 

23 

1.887 

535.2 

2.997 

544.7 

1.021 

523.1 

1.207 

525.9 

24 

1.988 

535.9 

3.119 

545.5 

1.054 

523.8 

1.251 

526.6 

25 

2.053 

536.6 

3.155 

546.3 

1.036 

524.5 

1.301 

527.4 

26 

2.122 

537.3 

3.250 

547.0 

1.112 

525.2 

1.311 

528.1 

27 

2.176 

538.0 

3.359 

547.8 

1.145 

525.8 

1.351 

528.8 

28 

2.207 

538.6 

3.455 

548.5 

1.159 

526.5 

1.415 

529.5 

29 

2.236 

539.2 

3.674 

549.2 

1.193 

527.1 

1.443 

530.1 

30 

2.255 

539.9 

3.812 

549.8 

1.217 

527.7 

1.468 

530.7 

31 

2.277 

540.5 

3.856 

550.5 

1.246 

528.3 

1.484 

531.4 

32 

2.327 

541.1 

3.726 

551.1- 

1.282 

528.9 

1.526 

532.0 

33 

2.384 

541.7 

3.886 

551.7 

1.313 

529.4 

1.578 

532.6 

34 

2.459 

542.3 

3.924 

552.3 

1.334 

530.0 

1.618 

533.2 

35 

2.537 

542.9 

4.009 

552.9 

1.367 

530.5 

1.670 

533.7 

36 

2.608 

543.5 

4.060 

553.5 

1.393 

531.1 

1.672 

534.3 

_ 37 _ 

2.688 

_544.0 

4.188 

554.1 

1.413 

531.6 

1.710 

534.9 

33 

2 . 754 

544.5 

4.222 

554.7 

1.447 

532.1 

1.749 

535.4 

39 

2.826 

545.0 

4.286 

555.2 

1.469 

532.6 

1.796 

535.9 

40 

2.875 

545.5 

4.292 

555.8 

1.489 

533.1 

1.788 

536.4 

41 

2.912 

546.0 

4.366 

556.4 

1.504 

533.6 

1.792 

536.9 

42 

2.940 

546.5 

4.400 

556.9 

1.52  5 

534.0 

1.827 

537.5 

43 

2.989 

546.9 

4.493 

557.4 

1.53  8 

534.5 

1.845 

538.0 

44 

3.014 

547.4 

4.537 

558.0 

1.542 

535.0 

1.846 

538.5 

45 

3.025 

547.9 

4.603 

558.5 

1.554 

535.4 

1.854 

539.0 

46 

3.032 

548.3 

4.628 

559.0 

1.564 

535.9 

1.827 

539.5 

47 

3.035 

548.8 

4.602 

559.5 

1.564 

536.4 

1.861 

540.0 

48 

3.048 

549.3 

4.706 

560.0 

1.564 

536.8 

1.876 

540.5 

49 

3.113 

549.7 

4.783 

560.5 

1.564 

537.3 

1.865 

541.0 

50 

3.114 

550.2 

4.896 

561.0 

1.549 

537.7 

1.856 

541.5 
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PERCENT 

RATE  OF 

TEMP. 

RATE  OF 

TEMP. 

RATE  OF 

TEMP. 

RATE  OF 

TEMP. 

WT.LOSS 

WT.LOSS 

DEG  C 

WT.LOSS 

DEG  C 

WT.LOSS 

DEG  C 

WT.LOSS 

DEG  C 

51 

3.128 

550.6 

4.917 

561.5 

1.542 

538.2 

1.869 

542.0 

52 

3.106 

551.1 

4790^3 

561.9 

1 .53  5 

538.7 

1.874 

542.4 

53 

3.108 

551.5 

4.951 

562.4 

1.521 

539.2 

1.858 

542.9 

54 

3.106 

552.0 

4.978 

562.9 

1.514 

539.6 

1.838 

543.4 

55 

3.125 

552.4 

5.055 

563.4 

1.517 

540.1 

1.832 

543.9 

56 

3.139 

552.9 

5.041 

563.8 

1.513 

540.6 

1.827 

544.4 

57 

3.146 

553.3 

5.047 

564.3 

1.514 

541. C 

1.898 

545.0 

58  " 

■  3.111 

553.8 

5.120 

564.8 

1.506 

541T5" 

1.92T^ 

5^.4 

59 

3.116 

554.2 

5.144 

565.2 

1.504 

542.0 

1.930 

545.9 

60 

3.100 

554.7 

5.196 

565.7 

1.495 

542.5 

1.890 

546.4 

61 

3.136 

555.  1 

5.185 

566.1 

1.480 

543.0 

1.810 

546.9 

62 

3.125 

555.6 

5.150 

566.6 

1.463 

543.4 

I  •  81 6 

547 .4 

63 

3.113 

556.0 

5.190 

567.0 

1.448 

543.9 

1.776 

547.9 

64 

3.086 

556.4 

5.179 

567.5 

1.434 

544.4 

1.75? 

548 . 4 

65 

3.052 

556.9 

5.215 

567.9 

1.425 

544.9 

1.732 

549.0 

66 

3.014 

557.4 

5.207 

568.4 

1.429 

545.4 

1.763 

549.5 

67 

2.980 

557.8 

5.138 

568.8 

1.440 

546.  C 

1.779 

550.0 

68 

2.948 

558.3 

5.221 

569.3 

1.437 

546.4 

1.767 

550.5 

69 

2.935 

558.8 

5.245 

569.7 

1.419 

546.9 

1.714 

551.0 

70 

2.920 

559.3 

5.270 

570.1 

1.390 

547.4 

1.705 

551.6 

71 

2.895 

559.8 

5.268 

570.6 

1.379 

547.9 

1.707 

552.1 

72 

2.862 

560.2 

5.175 

571.0 

1.357 

548.5 

1.716 

552.6 

73 

2.890 

560.7 

5.108 

571.5 

1.348 

549.0 

1.721 

553.2 

74 

2.880 

561.2 

5.065 

571.9 

1.361 

549.5 

1.696 

553.7 

75 

2.861 

561.7 

4.968 

572.4 

1.368 

550.1 

1.686 

554.2 

76 

2.846 

562.2 

4.933 

572.9 

1.368 

550.6 

1.683 

554.8 

77 

2.829 

562.7 

4.862 

573.4 

1.368 

551.1 

1.688 

555.3 

78 

2.779 

563.2 

4.914 

573.8 

1.368 

551.6 

1.743 

555.9 

79 

2.721 

563.7 

4.949 

574.3 

1.368 

552.1 

1.761 

556.4 

80 

2.686 

564.2 

4.931 

574.8 

1.368 

552.7 

1.739 

556.9 

81 

2.631 

564.7 

4.830 

575.2 

1.368 

553.2 

1.730 

557.4 

82 

2.577 

565.3 

4.666 

575.7 

1.368 

553.7 

1.680 

557.9 

83 

2.508 

565.8 

4.476 

576.3 

1.368 

554.2 

1.668 

558.5 

84 

2.514 

566 .4 

4.350 

576.8 

1.368 

554.7 

1.637 

559.0 

85 

2.475 

567.0 

4.204 

577.3 

1.354 

555.2 

1.627 

559.6 

86 

2.443 

567.5 

4.093 

577.9 

1.354 

555.8 

1.632 

560.1 

87 

2.411 

568.1 

3.996 

578.5 

1.347 

556.3 

1.614 

560.7 

88 

2.349 

568.7 

3.851 

579.1 

1.326 

556.8 

1.578 

561.2 

89 

2.278 

569.3 

3.734 

579.7 

1.303 

557.4 

1.539 

561.8 

90 

2.218 

569.9 

3.574 

580.4 

1.272 

557.9 

1.489 

562.4 

91 

2.145 

570.5 

3.442 

581.0 

1.235 

558.5 

1.420 

563.0 

92 

2.049 

571.3 

3.188 

581.8 

1.196 

559.1 

1.340 

563.7 

93 

2.009 

571.9 

3.016 

582.5 

1.158 

559.7 

1.255 

564.3 

94 

1.922 

572.6 

2.845 

583.4 

1.104 

560.3 

1.160 

565.1 

95 

1.818 

573.4 

2.652 

584.2 

1.028 

560.9 

1.059 

565.9 

96 

1.664 

574.1 

2.459 

585.2 

0.923 

561.6 

0.946 

566.8 

97 

1.454 

575.1 

2.180 

586.2 

0.794 

562.5 

0.817 

567.8 

98 

1.181 

576.1 

1.783 

587.5 

0.608 

563.5 

0.649 

569.0 

99 

0.827 

577.6 

1.115 

589.2 

0.392 

564.9 

0.429 

570.6 
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XXX 
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